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Research  initiated  under  contract  N00014-78-C-0889  is  most  conveniently 
summarized  in  the  form  of  four  appendices.  First,  in  Appendix  I,  the  five 
different  areas  of  research  on  water  waves  are  briefly  summarized.  In  two  of 
these  areas  (evolution  of  long  internal  waves  and  evolution  of  envelope  solitons) 
work  has  not  been  completed;  some  results  to  date  are  described.  In  the  remaining 
three  areas,  final  manuscripts  have  been  completed  and  are  presented  in  Appendices 
II,  III,  and  IV.  In  particular,  Appendix  II  is  a  copy  of  a  paper  entitled 
'Baroclinic  Tsunami  Generation"  to  appear  in  the  Journal  of  Physical  Oceanography , 
September  1980;  Appendix  III  is  a  copy  of  a  paper  entitled  "Long  Waves  Generated 
by  Complex  Bottom  Motions"  to  appear  in  Proceedings ,  17th  Conference  on  Coastal 
Engineering;  and  Appendix  IV  is  a  copy  of  lecture  notes  on  small-scale  ocean 
waves  presented  at  the  summer  (1980)  course  "Topics  In  Ocean  Physics"  held  at 
the  International  School  of  Physics  "Enrico  Fermi",  Varenna,  Italy. 


Appendix  I 
Research  Summary 


Research  Summary 


Research  conducted  during  the  past  year  under  Contract  No.  N00014-78-C-0889 
may  be  conveniently  divided  into  four  categories. 

I.  Baroclinic  Tsunami  Generation.  Utilizing  experimental  data  obtained 
previously,  a  study  was  completed  on  the  internal  waves  generated  in  a  stratified 
ocean  by  a  vertical  motion  of  the  sea  floor.  Details  of  this  study  are  presented 
in  the  enclosed  manuscript  (see  Appendix  II)  which  will  appear  in  the  Journal  of 
Physical  Oceanography. 

II.  Evolution  of  Long  Internal  Waves.  Subsequent  to  the  generation  processes 
mentioned  above,  the  evolution  of  long  internal  waves  on  a  thin  pycnocline  has 

been  examined.  Dr.  Harvey  Segur  of  Aeronautical  Research  Associates  of  Princeton 
and  I  have  analyzed  the  experimental  data  and  prepared  a  draft  manuscript  which 
examines  several  model  equations  for  internal  wave  evolution  in  the  context  of 
experimental  data.  This  work  is  ongoing;  a  copy  of  our  final  manuscript  will 
be  forwarded  as  soon  as  possible. 

III.  Long  Waves  Generated  by  a  Complex  Bottom  Motion.  As  a  final  statement 
on  barotropic  water-wave  generation  by  motions  of  the  sea  floor,  Professor  Frederic 
Raichlen  (of  Cal  Tech)  and  I  completed  a  manuscript  (enclosed)  which  illustrates 

a  technique  for  computing  wave  structure  when  the  time-displacement  history  of 
the  sea  floor  is  very  complicated.  The  analytical  mode  was  verified  by  experiments. 
These  results  were  presented  at  the  17th  International  Conference  on  Coastal 
Engineering  in  Sydney,  Australia  (March  1980) ,  and  the  manuscript  (see  Appendix  III) 
will  appear  in  the  meeting  proceedings. 


IV.  Lecture  Notes  on  Small-Scale  Ocean  Waves.  At  the  invitation 


of  the  Italian  Society  of  Physcis,  I  presented  a  series  of  lectures  at  their 
summer  course  on  "Topics  in  Ocean  Physics"  held  during  7  July  -  19  July  1980  at 
the  International  School  of  Physics  "Enrico  Fermi",  in  Varenna,  Italy.  These 
lectures  focused  on  the  nonlinear  aspects  of  water  waves  including  solitons  and 
resonant-interaction  forcing  of  edge  waves.  The  collection  of  lecture  notes 
from  the  ten  invited  speakers  will  be  published  by  the  North  Holland  Publishing 
Company.  A  copy  of  my  notes  is  presented  in  Appendix  IV. 

V.  Evolution  of  Short-Wave  (Envelope)  Solitons.  The  first  goal  of  the 
experimental  research  program  on  short-wave  evolution  was  the  direct  generation 
of  an  envelope  soliton.  The  best  results  obtained  to  date  for  a  wave  packet 
with  a  carrier  wave  period  T  =  1  sec  and  water  depth  h  =  1  m  are  shown  in  figures 
1  abed.  Results  at  four  positions  along  the  (Berkeley)  wave  tank  are  presented 
with  theoretical  soliton  profiles  based  on  the  local  amplitude  of  the  wave  packet. 
The  initial  data  measured  at  x  =  6  m  from  the  wave  paddfe  agree  well  with  the 
theoretical  soliton  profile.  During  subsequent  propagation  the  packer  retains 
its  compactness  and  travels  with  the  (predicted)  linear  group  speed  of  the  carrier 
wave.  However,  the  measured  packet  deviates  slightly  from  the  predicted  soliton 
profile;  this  behavior  is  especially  obvious  at  the  last  station  (x  =  45  m)  of 
measurement.  The  front-to-back  asymmetry  emerging  at  x  =  45  m  is  similar  to  that 
observed  by  previous  researchers  on  deep  water  wave  packets.  In  fact,  the 
symmetrical  measurements  prior  to  the  last  measurement  station  are  quite  exceptional 
and  have  not  been  reported  previously  in  the  literature.  (I  might  add  that  the 
first  three  wave  traces  of  figure  1  are  almost  identical  to  those  generated  in 
the  Florida  facility;  only  the  added  length  of  the  Berkeley  wave  tank  permitted 
the  distortion  to  be  clearly  observed.) 
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Below  each  wave  trace  in  figure  1  I  have  also  presented  a  histogram 
of  the  node-to-node  wave  periods  through  the  wave  packet.  This  rather 
crude  attempt  at  complex  demodulation  provides  some  insight  into  the  wave 
packet  evolution.  At  the  first  two  stations,  local  wave  periods  differ 
little  from  the  programmed  period  of  1  sec.  However,  at  the  last  two 
stations  the  leading  waves  are  clearly  developing  periods  up  to  15%  longer 
than  the  1  sec  period  of  waves  near  the  center  of  the  packet.  The  origin 
of  the  longer  wave  content  is  not  clear  although  its  emergence  at  the 
front  of  the  packet  is  probably  responsible  for  the  changes  between  meas¬ 
ured  and  predicted  packet  piofiles.  The  most  plausible  explanation  for 
the  unwanted  wave  content  of  the  packet  is  that  it  is  a  consequence  of 
imperfect  generation  by  the  wave  maker.  It  is  an  unfortunate  (theoretical 
and  observed)  property  of  short  wave  solitons  that  these  nonlinear  waves 
propagate  with  the  linear  group  speeds  of  their  carrier  waves.  Hence, 
unlike  their  long-wave  counterparts,  short-wave  solitons  cannot  separate 
from  any  background  "noise"  with  the  same  frequency.  (The  background 
noise  corresponds  to  the  radiation  or  non-sol iton  component  of  the  solution 
to  the  nonlinear  Schrodinger  equation.)  The  mergence  and  slow  migration 
(at  linear  speeds)  of  the  longer  wave  components  in  figure  1  may  he  a 
i  ,ui i fest at i on  of  this  phenomenon.  In  this  regard  it  is  important  to  note 
that  the  radiation  wave  components  remain  dispersive  even  at  nonlinear 
order  and  their  amplitudes  should  decay  in  time  t  with  an  inviscid  decay 

-h 

rate  of  t  .  Since  the  inviscid  sol iton  does  not  decay,  when  noise  is 
present  one  expects  the  measured  data  to  agree  hotter  with  the  theoretical 
profile  with  distance  (or  time)  down  the  tank.  fxactly  the  opposite 


8 


bebavi  or  is  observed  in  figure  1;  however,  it  mist  be  remembered  that  glow¬ 
ing  distortions  might  still  occur  until  the  radiation  attains  its  asymp¬ 
totic  character. 

Although  inviseid  solitons  do  not  decay  with  time,  water-wave  solitens 
are  subject  to  viscous  attenuation  of  their  amplitude.  Two  models  for  the 
viscous  damping  of  short-wave  solitons  have  been  developed  independently  by 
Dr.  Harvey  Segur  and  myself  (with  considerable  assistance  from  N.  K.  Lin, 
a  graduate  student  at  Berkeley).  Both  models  predict  similar  results  even 
though  the  approaches  to  model  development  were  quite  different.  Experi¬ 
mental  verification  of  the  models  must  be  considered  limited  for  large  pro¬ 
pagation  distances  due  to  the  distortion  discussed  previously;  short  distance 
evolution  according  to  these  models  is  verified  by  the  measured  data.  As 
in  all  \ i  •"  mis  damping  models,  contamination  of  the  water  surface  enters  as 
an  e  pi  l  ie.  a  1  loefficient  ranging  from  zero  for  no  contamination  to  unity 
for  a  fully  contaminated  (and  immobile)  surface.  Theoretical  results  for 
these  two  limiting  cases  are  shown  in  figure  2  with  experimental  data  for 
the  wave  evolution  in  figure  1.  Measurements  suggest  that  some  surface 
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c  on  t  a  i  n  a  t  i  on  is  present  in  the  Berkeley  data  which  is  expected  since  water 
:iins  in  the  tank  continuously.  Theoretical  and  experimental  results 
for  the  ‘-.esc  initial  data  observed  in  the  Florida  facility  are  also  shown 
in  figure  2.  In  Florida,  well  water  was  pumped  to  the  wave  tank  immediately 
prior  to  the  experiments;  lance,  surface  contamination  should  he  less;  this 


expectation  is  verified  by  the  data.  (Ihc  narrower  tank  in  Florida  leads 
to  greater  wave  attenuation  with  distance  of  propagation.) 

Tt  should  also  be  noted  that  the  viscous  damping  models  for  envelope 
solitons  indicate  that  they  decay  twice  as  fast  as  a  linear^  wave  train  of 
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the  same  frequency.  Unfortunately,  water  wave  trains  of  sufficiently  small 
amplitude  were  not  measured  to  test  this  prediction.  However,  measurements 
of  wave  trains  with  the  same  degree  of  nonlinearity  as  the  comparison  soli- 
ton  (prior  to  the  visual  onset  of  Benjamin-Feir  instability)  indicate  that 
the  solitons  decay  only  one-half  as  rapidly  as  the  uniform  wave  train. 

The  evolution  of  arbitrary  packets  of  initial  data  was  also  investi¬ 
gated.  An  example  of  results  is  shown  in  figures  3a-b  at  two  measurement 
stations;  histograms  of  wave  periods  are  shown  in  figures  4a-b  for  each 
measurement.  The  dominant  wave  packet  in  the  downstream  measurement  of 
figure  3b  clearly  has  evolved  a  soliton  profile.  The  remaining  wave  packets 
do  not  have  soliton  profiles  and  represent  the  radiation  component  in  the 
initial  data.  The  period  histogram  of  figure  4b  clearly  shows  the  linear 
dispersive  nature  of  the  constituent  wave  components  which  are  rank  ordered 
by  their  group  speed.  (Also  see  Appendix  IV) . 

Experiments  on  the  stability  of  2-D  envelope  solitons  to  3-D  perturba¬ 
tions  have  proven  to  be  inconclusive.  Very  small  3-D  effects  were  measured 
in  some  experiments  but  without  a  consistency  or  magnitude  sufficient  to 
infer  definitive  behavior  patterns.  The  experiments  did  highlight  a  defi¬ 
ciency  in  the  analytical  studies  on  this  aspect  of  envelope  soliton.  All 
explicit  criteria  which  provide  definitive  predictions  on  the  tank  (or 
crest)  widths  necessary  for  instability  assume  infinite  water  depths.  How¬ 
ever,  the  nonlinear  description  of  these  waves  contains  terms  which  decay 
algebraically--not  exponentially--with  depth.  Hence,  it  is  dangerous  to 
apply  these  criteria  to  finite  depth  experiments,  even  for  kh  ~  10. 
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Finally,  based  on  the  data  presented  here  as  well  as  many  other  experi¬ 
ments  on  varying  aspects  of  short-wave  solitons,  it  is  clear  to  me  that  the 
reality  of  these  systems  is  more  difficult  to  disentangle  than  that  of  their 
long-wave  counterparts.  However,  even  with  the  confusion  arising  from 
microscopic  observation  of  these  waves,  it  is  also  apparent  that  the  macro¬ 
scopic  structure  predicted  by  the  nonlinear  Schrodinger  equation  is  present. 
All  of  these  results  emphasize  the  need  for  more  analytical  study  coupled 
with  delicate  experiments  capable  of  yielding  quantitative  measurement  of 
third  order  (and  higher)  parameters. 
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Baroclinic  Tsunami  Generation 

Joseph  L.  Hammack 

Department  of  Civil  Engineering,  University  of  California.  Berkeley  94720 
(Manuscript  received  13  August  1979.  in  final  form  2  May  19K0i 

ABSTRACT 

An  analytical  and  experimental  study  of  the  baroclinic  waves  generated  bt  a  monopole  dislocation 
of  the  sea  floor  is  presented.  Analytical  results  are  based  on  a  two-dimensional  and  linearized  description 
of  motion  using  a  two-layer  approximation  for  density  variation;  experiments  utilize  a  stratification 
with  finite  (nonzero)  pycnocline  thickness.  Scaling  parameters  which  characterize  the  generation  process 
and  the  potential  role  of  nonlinear  effects  are  discussed.  It  is  shown  that  the  barotropic  modes  are 
not  affected  by  the  small  differences  in  potential  density  typical  of  ocean  stratifications  and  all  previous 
results  for  these  waves  are  applicable.  The  two-layer  approximation  is  found  to  provide  an  accurate  repre¬ 
sentation  of  the  (long)  baroclinic  waves  typical  of  tsunamis.  Like  the  barotropic  response,  baroclinic 
generation  is  impulsive  and  linear  resulting  in  wave  amplitudes  proportional  to  the  vertical  offset  of  the 
sea  floor.  Near  the  generation  region  barotropic  waves  have  amplitudes  of  one-half  the  sea  floor  displace¬ 
ment  while  the  baroclinic  waves  are  attenuated  further  by  the  ratio  hjh,  where  h  is  the  total  fluid  depth  and 
h,  the  upper  layer  thickness.  Although  Coriolis  effects  are  not  included  in  either  the  analytical  or  experi¬ 
mental  models,  these  effects  may  often  be  significant  for  baroclinic  waves.  In  general,  the  potential  role 
of  Coriolis  forces  is  both  earthquake  and  site  specific  Regardless,  the  analysis  herein  remains  applicable 
for  times  smaller  than  the  local  inertial  period 


1.  Introduction  and  major  conclusions 

Tsunami  research  to  date  has  generally  focused  on 
the  waves  at  the  ocean's  surface,  ignoring  vertical 
variations  in  potential  density  which  typically  exist. 
Submarine  earthquakes  which  produce  significant 
surface  (barotropic)  waves  will  also  generate  in¬ 
ternal  (baroclinic)  tsunamis  in  stratified  regions. 
These  baroclinic  modes  were  examined  first  by 
Cherkesov  (1968)  who  presented  a  theoretical  in¬ 
vestigation  of  waves  generated  by  a  motion  of  the  sea 
floor  in  an  ocean  with  an  exponential  increase  in  po¬ 
tential  density  with  depth.  Subsequently.  Fedo- 
senko  and  Cherkesov  (1968)  presented  similar  re¬ 
sults  for  a  two-laver  stratification  which  is  more 
representative  of  actual  ocean  conditions.  Solutions 
were  based  on  a  linearized  description  of  motion  for 
both  two-dimensional  and  three-dimensional  (axi- 
symmetric)  source  regions;  Coriolis  forces  were 
neglected.  In  both  presentations,  the  integral  solu¬ 
tions  were  evaluated  only  for  the  far-field  at  large 
distances  from  the  source  region.  These  results 
showed  that  a  displacement  of  the  sea  floor  pro¬ 
duced  (fast)  barotropjc  waves  and  (slow)  baro¬ 
clinic  waves  which  separated  with  time.  Each  wave 
system  induced  a  coupled  motion  at  the  opposite 
interface.  It  was  further  claimed  that  the  asymp¬ 
totic  amplitudes  of  the  internal  waves  could  exceed 
those  at  the  free  surface  by  a  factor  of  2?  for 
ocean  conditions.  Neither  the  relevance  of  asymp¬ 
totic  linear  analysis  to  actual  tsunamis  nor  any 
limitation  of  the  results  due  to  omission  of  Coriolis 
forces  were  discussed  by  the  authors.  [Hammack 
and  Segur  (1978)  have  shown  that  linear  asymp¬ 
totic  solutions  never  become  applicable  for  baro¬ 
tropic  tsunamis.) 
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The  objective  of  the  present  paper  is  to  examine 
the  near-field  behavior  of  baroclinic  waves  gener¬ 
ated  by  vertical  motions  of  the  sea  floor.  Both 
analytical  and  experimental  results  are  presented 
for  a  simple  model  of  generation  consisting  of  a 
block  upthrust  (or- downthrow)  according  to  a  pre¬ 
scribed  time-displacement  history.  Analytical  re¬ 
sults  are  based  on  a  two-layer  stratification  of 
potential  density  while  the  experiments  utilize  a 
stratification  with  finite  (nonzero)  pycnocline 
thickness. 

The  organization  of  this  paper  and  some  of  the 
major  results  are  as  follows.  In  Section  2  a  brief 
description  of  tsunamigenic  earthquakes  is  pre¬ 
sented  and  the  model  to  be  examined  herein  is 
introduced.  Nondimensional  parameters  which 
characterize  the  generation  process  are  discussed 
in  Section  3  where  it  is  shown  that  the  generation 
process  of  actual  tsunamis  is  impulsive  and  linear 
for  both  the  barolropic  and  baroclinic  response.  In 
Section  4  the  linear  solution  for  the  baroclinic 
model  is  developed  and  other  well-known  properties 
of  long  internal  waves  at  the  interlace  of  two 
liquids  are  reviewed.  In  the  Boussinesq  limit  of 
small  density  differences,  it  is  shown  that  the  baro- 
tropic  response  is  not  affected  by  the  stratification. 
In  addition,  both  barotropic  and  baroclinic  tsunamis 
have  initial  shapes  similar  to  the  permanent  de¬ 
formation  of  the  sea  floor.  The  maximum  amplitude 
of  the  surface  wave  is  one-half  the  vertical  dis¬ 
placement  of  the  sea  floor  while  the  maximum 
amplitude  of  the  internal  wave  is  additionally 
attenuated  by  the  ratio  of  upper  layer  depth  to  total 
depth  Experimental  equipment  and  procedures  are 
described  in  Section  5  and  theoretical  results  are 
compared  with  the  experimental  measurements  in 
Section  6.  Results  indicate  that  the  maximum 
amplitudes  of  the  barotropic  and  baroclinic  waves 
near  a  source  region  may  be  represented  by  a  single 
(and  simple)  functional  relation  over  the  full  range  of 
their  respective  generation  time  scales  (see  Fig.  10). 
The  potential  role  of  Coriolis  effects  for  baroclinic 
waves  is  discussed  in  Section  7,  where  it  is  dem- 
onstiated  that  time  scales  of  baroclinic  modes  can 
a  1  nui)  often  approach  the  inertial  period. 


2.  Tsunamigenic  earthquakes  and  the  proposed 
model 

Most  large  tsunamis  are  associated  with  sub- 
oceanic  earthquakes,  centered  near  plate  margins, 
during  which  significant  vertical  offset  of  the  sea 
floor  occurs.  In  order  to  determine  the  waves 
generated  by  mechanisms  of  this  type,  a  knowledge 
of  the  water  -depths,  stratification  and  latitude  at 
the  source  site  as  well  as  the  spatial  and  tem¬ 
poral  distributions  of  the  sea  floor  dislocation  is 
required.  Since  the  objective  herein  is  to  delineate 
fundamental  relations  between  the  generation 
mechanism  and  wave  structure,  we  seek  a  tractable 
model  which  embodies  essential  source  parameters 
in  a  simple  way.  With  this  goal  in  mind  we  will 
assume  a  uniform  depth  It  in  the  generation  re¬ 
gion  and  adopt  a  two-layer  model  to  represent 
the  ambient  stratification  in  potential  density.  In 
addition  we  will  neglect  Coriolis  forces  (and  their 
latitude  dependence);  limitations  of  the  results  due 
to  this  omission  will  be  discussed  a  posteriori. 
Before  introducing  further  simplifications,  a  brief 
description  of  features  common  to  tsunamigenic 
earthquakes  and  their  (dimensional)  scales  is 
presented. 

It  is  both  convenient  and  appropriate  to  define  the 
source  size  of  the  tsunami  (and  earthquake)  as 
the  area  affected  by  aftershocks  immediately 
following  the  generative  earthquake.  This  area  is 
typically  elongated  and  elliptical  in  shape  with  the 
ruptured  portion  of  the  fault  joining  the  foci  of  the 
ellipse.  Early  efforts  to  correlate  the  source  area 
S  with  Richter  magnitude  1  Ms  proved  unsatisfactory 
due  to  the  marked  scatter  in  the  data.  Kanamori 
and  Anderson  (1975)  inter  alios  have  shown  that 
definitive  correlations  exist  between  S  (km2)  and 
the  seismic  moment  Mn  (dyn  cm).  Representative 
values  based  on  their  results  for  interplate  earth¬ 
quakes  are  shown  in  Table  1.  In  order  to  provide 
perspective,  it  is  still  desirable  to  relate  the  seismic 
moment  to  Richter  magnitude.  As  expected,  results 
here  are  less  definitive  and  no  single  empirical 
relation  exists;  however,  Kanamori  and  Anderson 
(1975)  do  present  correlations  for  three  classes  of 
earthquake  magnitudes.  Approximate  values  of  M , 
based  on  their  results  are  shown  in  Table  I.  As  a 
first  approximation,  we  may  assume  that  the  source 
areas  of  Table  1  represent  regions  of  monopole 
vertical  offset  (upthrust  of  downthrow).  With  this 
simplification  the  source  size  can  be  characterized 
by  a  single  length  dimension  L  ~  51'2;  representative 
values  are  shown  in  Table  I.  If  we  further  assume 
the  static  and  seismic  moments  are  equal,  mag¬ 
nitudes  of  the  average  (and  permanent)  vertical 
offset  £„  can  be  found  by 
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!C,,|  =  A iJyS, 

where  n  is  the  shear  modulus  (rigidity)  of  the  frac¬ 
tured  medium.  Estimates  ot  |£„|  with  p.  =  5  x  10" 
dynes  cm  2  are  presented  in  Table  1.  The  temporal 
history  of  the  sea  floor  dislocation  is  characterized 
by  a  minimum  of  two  time  scales:  a  time  tc  cor¬ 
responding  to  the  time  required  for  the  vertical 
offset  to  occur  at  a  specific  spatial  location  and  a 
time  t„  representing  the  period  required  for  the 
rupture  to  propagate  along  the  fault.  Empirical 
results  of  Kanamori  and  Anderson  (1975)  suggests 
that  tH  ~  10/r  for  tsunamigenic  earthquakes  (Ms 
s  6.5)  with  estimates  of  the  rupture  time  ranging 
from  tR  -  10  s  for  moderate  earthquakes  to  tR 
~  100  s  for  earthquakes  with  fault  lengths  exceeding 
100  km  (based  on  rupture  speeds  VR  ~  LltR  ~  2-3 
km  s'1).  As  suggested  previously,  vertical  offset  is 
primarily  responsible  for  tsunami  generation;  hence, 
the  time  tr  is  most  relevant  for  subsequent 
scaling. 

In  light  of  the  above  description,  we  will  adopt 
the  following  simplifications  for  the  generation 
model.  First,  the  smallness  of  the  linear  source 
dimensions  of  Table  1  relative  to  the  earth  s 
radius  permits  a  plane-earth  approximation.  Second, 
for  monopole  and  uniform  vertical  offset,  the 
elongation  of  the  source  dimension  suggests  that 
initial  propagation  of  waves  near  the  center  of  the 
source  region  is  one-dimensional  and  normal  to  the 
fault.  Hence,  we  arrive  at  the  simple  generation 
model  shown  in^  Fig.  I.  The  actual  sea  floor  dis¬ 
placement  used  in  computations  is  given  by 

£(.r ./)  =  ^„(e  -  1  )H{b1  -  *2),  (1) 

where  H  represents  the  Heavyside  stepfunction.  In 
the  dislocation  of  ( I ),  we  have  taken  2b  ~  L  and,  for 
experimental  convenience,  we  will  take 


a 


3.  Scaling  the  generation  process  and  a  comment  on 
nonlinearity 

In  order  to  provide  a  rational  basis  for  choosing 
computational  examples,  nondimensional  param¬ 
eters  which  characterize  generation  and  typical 
magnitudes  for  tsunamigenic  earthquakes  are  re¬ 
quired.  With  this  objective  in  mind,  a  brief  review 
of  results  by  Hammack  (1973)  for  generation  of 
barotropic  waves  is  useful.  Three  nondimensional 
parameters  are  required  to  characterize  barotropic 
generation  for  simple  models  similar  to  that  dis¬ 
cussed  in  Section  2:  an  amplitude  scale.  |(„| Ih:  a 
size  scale,  bllr.  and  r  =  trC,/b  termed  a  time- 
size  ratio  where  C,  =  (gh)1 2  is  the  speed  of  long 
barotropic  modes.  [Alternatively,  r  may  be  con¬ 
sidered  an  inverse  Froude  number  analogous  to  that 
used  by  Noda  ( 1970)  for  landslide-generated  waves.) 
The  generation  process  may  be  classified  according 
to  the  time-size  ratio  as  impulsive  for  t4  I,  creep¬ 
ing  for  t  >  1  and  transitional  for  r  =  1.  The  specific 
range  of  t  for  each  classification  is  only  a  weak 
function  of  the  size  scale  with  the  transitional  regime 
essentially  disappearing  for  large  size  scales  (b/h 
s>  1).  The  amplitudes  of  waves  propagating  out  of 
the  generation  region  for  an  impulsive  motion  are 
proportional  to  the  bed  displacement  ( rj  cc  („) 
where  the  constant  of  proportionality  depends 
only  on  the  size  scale.  For  b/h  >  I,  this  constant 
attains  its  maximum  value  of  one-half  and  the 
tsunami  near  the  source  region  is  similar  in  shape  to 
the  permanent  deformation  of  the  sea  floor.  For 
creeping  bed  motions,  wave  amplitudes  are  reduced 
in  inverse  proportion  to  the  magnitude  of  the  time- 
size  ratio  (17  °c  £0 It).  Nonlinear  effects  become 
significant  when  generation  is  impulsive  and 
3  0.2.  Assuming  depths  in  the  generation  region 
of  h  =  200  m  on  the  shelf  and  h  *=  4000  m  in  deep 
water,  size  scales  for  tsunamigenic  earthquakes  are 
typically  large  (based  on  the  lengths  L  ~  b  given 
in  Table  1);  hence,  energy  is  concentrated  at  long 
wavelengths.  In  addition,  amplitude  scales  are  small 
based  on  values  of  |£„|  given  in  Table  I.  Finally, 
using  characteristic  times  up  to  t,.  —  10  s.  the  baro¬ 
tropic  response  is  typically  impulsive  (-  <S  1 ).  Based 
on  these  characteristic  values,  the  barotropic  re¬ 
sponse  is  expected  to  be  linear  and  the  wave  struc¬ 
ture  near  the  source  region  should  resemble  the 
permanent  deformation  of  the  sea  floor  with  a 
50%  attenuation  in  amplitude.  It  is  noteworthy,  and 
fortunate,  that  the  barotropic  response  is  repre¬ 
sented  by  the  simplest  possible  limiting  conditions 
of  the  generation  parameters. 


Physically.  tc  represents  the  time  required  for  two- 
thuds  of  the  vertical  offset  to  occur. 
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An  examination  of  the  independent  variables 
( p1,p2,/i,,A2,A,r,.,£u,g)  for  generation  in  a  two-layer 
ocean  (see  Fig.  I)  indicates  that  five  nondimen- 
sional  parameters  are  required  to  characterize 
generation.  A  possible  parameter  choice  for  small 
density  differences  is 

blh 

7,  =  t,.CJb  =  Inf’ll)' Vb  ,  (3) 

r-j  =  irCtlh  =  IAeghJiJh)'-/b  I 

(/t, //I,)'1  j 

where  h  -  /i ,  +  /i->  and  e  =  (p2  -  Pi)/p->  1.  Of 

these  five  parameters,  the  first  four  preserve  the 
scaling  characteristics  of  the  barotropic  mode  and 
indicate  the  expectation  of  a  similar  role  for  a  time- 
size  ratio  (r-i )  based  on  the  internal  (long)  wave 
speed  C>  (see  Section  4).  Notice  that  the  two 
time-size  ratios  satisfy  t-Jt ,  =  0(e' 2);  hence,  bed 
motions  which  are  classified  creeping  for  the  surface 
response  (t,  ^  1)  may  be  impulsive  for  the  inter¬ 
facial  response  (r2  1).  Since  r,  <  1  for  prototype 

tsunamis,  the  baroclinic  response  is  always  impulsive. 

The  fifth  parameter  |(/i,//i2)  -  1)  is  geometric  in 
nature;  the  various  roles  of  this  parameter  and  its 
related  forms  A ,//r  or  /;.,//(  will  appear  in  Sections 
4  and  6.  Its  specific  role  as  an  indicator  of  the 
strength  of  nonlinear  effects  during  generation  can 
be  established  as  follows.  In  a  two-layer  system, 
nonlinear  effects  may  occur  for  both  the  surface 
and  interfacial  motion.  For  the  free  motion  of  long 
internal  waves  in  a  two-layer  system  with  small 
density  differences  and  uniform  depth,  it  can  be 
shown  (see  Keulegan,  1953)  that  the  appropriate 
measure  of  small  but  finite  nonlinearity  is  |t)2(/i, 
-  A2)//i,A2 |.  It  will  be  demonstrated  in  Section  4 
that  t),  cc  £0(/i,//i)  for  impulsive  generation;  hence, 
nonlinearity  is  measured  by  the  parameter  |(£0/ 
li)Uitili->  -  1)|.  Since  hjh>  =£  l  for  typical  ocean 
stratifications  and  (£0//i)  is  small,  both  barotropic 
and  baroclinic  free  motions  are  linear  initially. 

4.  Theoretical  analysis 

We  consider  a  two-dimensional  ocean  as  shown 
in  Fig.  I  consisting  of  two  fluid  domains  Dx  and 
D,  bounded  above  by  a  free  surface  S,,  below  by  a 
solid  boundary  5*  and  contacting  at  a  liquid-liquid 
interface  5,.  Both  fluid  regions  are  unbounded  in 
the  direction  a  of  wave  propagation.  The  upper  and 
lower  fluids  have  densities  p,  and  p2,  respectively, 
and  are  gravitationally  stable  (p,  <  p2).  Initially,  the 
two  fluids  are  at  rest  with  uniform  depths  li,  and 
A,.  Subsequently  (/  >0),  the  solid  boundary  S,, 
is  permitted  to  move  in  a  ptescribed  manner  given 
by  y  =  -bn  +  ((a,/),  where  £(*,/)  is  given  by  (1). 
The  resulting  motion  of  the  free  surface  ry , ( .v . / )  and 
interface  T?2(.r,t)  are  required.  Since  both  motions 
are  probably  linear  for  prototype  tsunamis,  a  first- 
order  (linear)  description  of  motion  is  employed 
from  the  outset. 
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Assuming  the  fluids  in  each  domain  to  be  incom¬ 
pressible  and  their  motion  irrotational,  velocity- 
potentials  4>,(x .y  j)  and  <j>.,(x,y,t)  exist  for  Z),  and 
D...  respectively.  Conservation  of  mass  requires 

r-<i>,  =0  in  D,,  (4a) 

V =  0  in  (4b) 

where  V-  is  the  Laplacian  operator.  Assuming  the 
fluids  to  be  inviscid,  the  kinematic  boundary  condi¬ 
tions  at  each  surface  are 

</>,.  =  T)ti  on  y  =  A,.  (5a) 

=  V>,  on  y  =  0.  (5b) 

By  further  neglecting  surface  energy  at  Sf,  the  condi¬ 
tion  for  constancy  of  pressure  across  and  along  the 
free  surface  is 

4>i,  +  8Vt  =  °  on  y  =  /i , .  (6a) 

(The  hydrostatic  pressure  p,g/t,  has  been  adsorbed 
into  <bi.)  Continuity  of  pressure  across  5,  requires 

+  £92  =  —  (<£,,  +  gr)2)  on  y  =  0.  (6b) 

P  2 

The  kinematic  condition  at  the  solid  boundary  is 

4>i,  =  £<  on  y  =  — /i2,  (7) 

where  it  is  assumed  that  motion  of  the  solid  boundary 
is  not  affected  by  fluid  motion. 

Combining  the  kinematic  and  dynamic  boundary 
conditions  at  Sf  and  St  yields 

<*i„  +  8<t>i,  =  0  on  y  =  A , ,  (8) 

^2,,  +  =  —  (d>,„  +  g<t>,  )  on  v  =  0.  (9) 

Pi 

Using  the  Laplace  transform  in  t  and  the  Fourier 
transform  in  v.  Eqs.  (4),  (8),  (9)  and  (7)  become 


^i„  _  k-4>,  =  0  in  £>,,  (10a) 

-  k'4>i  =  o  in  Z?2,  ( 10b) 

s‘<t> i  +  =  0  on  v  =  A,,  (10c) 

.t2<f>,  +  £</>,„  =  —  (s'-'<j>,  +  gi>,  )  on  v  =  0,  (lOd) 

Pi 

4>i,  =  si  on  y  =  -h2,  ( lOei 


where  the  tilde  of  a  function /(v.r)  indica'e* 
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/( A ,.v )  =  I  dx  (  e>kJ’e~’f(x,t)dt.  (11) 

The  function/U./)  is  recovered  from  ( 1 1)  by  invert¬ 
ing  the  Laplace  and  Fourier  transforms,  i.e.. 

1  P  n 

./(*./>  =  —  dk 
2 7T  L* 

lim  —  f  '  *-kT"s'j\k,s)ds  ,  (12) 

'  2m 

using  the  complex  inversion  integral  for  the  Laplace 
transform.  The  transformed  surface  and  interfacial 
waves  are 


V\  = 


on  v  =  //,. 


(13a) 


Vi  =  —  l(Pi/p*)«^i  -  </>*!  on  v  =  0.  (13b) 

«.t? 

Solution  of  (10a,b)  yields 

<i>,(A.v,.v)  =  A,  coshAy  +  A2  sinhAv,  (14a) 

<h,(k,y,s)  =  A-,  coshAy  4-  A,  sinhAv,  (14b) 

and  the  coefficients  A,  =  A,(A.s)  with  i  =  1,4  are 
found  from  Eqs.  OOc)-(lOe)  to  be 


A  i  = 


"  O 


A,  -  A, 
A  i 


Dk  coshA/i  , 

-v:,A 


•v-7-,  +  ,</A  1 , 


(IA  coshA/n, 


I  vs  +  cA7, 


(/A  coshA/i, 


.v‘pr7,  +  s-gk  +  e(  jt?A_  )27~ , ) . 


(15a) 


( 1 5b) 

(15c) 


with 


8  =  s4(  I  +  —  7,7, )  +  v2(7,  +  7,)i>A 

’  />2  ' 

+  t(gA)-r,T,,  ( 16a) 

T,  =  tanhAA,,  i  =  1.2.  (16b) 

Finally,  the  transformed  free  surface  and  internal 
waves  are  found  from  (13)  to  be 


a4 

(17a) 

tl  coshA/i  i  coshA/z2 

- -  IP  oWI. 

0  coshkh-z 

(17b) 

The  transform  of  the  bed  motion  given  by  (1)  is 


aA-.v)  =  2£0 


sinAA> 


a 


k  Li(s  +  a) 


(18) 


Substituting  (18)  into  (17),  performing  the  integra¬ 
tion  around  the  Bromwich  contour,  taking  only  the 
real  part  and  noting  that  the  integrands  are  even 
functions  of  A,  the  free  surface  and  interfacial 
waves  are  found  to  be 


T)  ,(*./) 


T  /•(*./) 


r  B,Bi 
coshA/i,  coshA/t., 

x  [e-“'-  B.,B,  +  B:rB6]dk,  (19a) 


B.B , 


coshA/to 


i  + 


5.5  4 


with 


B  i 

«2 

S, 

5.V 


</A. 


(19b) 


-2/,  sinA/i  cosA.v 
;rA[l  +  <p,/p2)7,72] 


(20a) 


«■* 

(or  +  ci>t2)(or  +  (a,3) 
ai,3(  or  +  ta22) 
of(ai|2  -  o>>3) 

1D|  . 

COSoj,/  +  —  Sinca,/, 
« 

ta/fu3  +  a>,2) 
or2(<o,-’  +  <a22) 

cosw2/  +■  —  sima./ , 
a 


(20b) 

(20c) 

(20d) 

(20e) 

(200 


and  u>,-  U  =  I,  2)  are  the  roots  of  the  dispersion 
relation 

(1  +  p,/p,7 ,T2)oj4  -  gk(T,  +  T2)oj- 


given  by 


+  e{gktt,T,  =0.  (21) 


!  (7-,  +  T,)  ±  |(T,  +  72)2 

,  j  -  4e(l  +  p,/p27,72)7,72]''2 

<uT.2  =  gA  j  -  .  (22) 

2(1  +  p,/p27,72) 

(Note  that  io  =  -is  in  the  above  analysis.)  Eq.  (2 1 )  is 
well  known  for  two-layer  systems  with  a,,2  corre 
sponding  to  the  barotropic  modes  and  o>22  corre 
sponding  to  the  baroclinic  modes  (see  e  g  I  n"4 
1932  p  37?) 
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Although  the  Fourier  integral  solutions  given  by 
(19)  are  algebraically  complex  and  the, wave  struc¬ 
ture  is  difficult  to  see,  it  is  easily  shown  that  these 
solutions  contain  more  familiar  results.  For  ex¬ 
ample,  in  the  Boussinesq  limit  of  small  density 


differences  <e 

— *  0)  we  find 

1  ‘  cos A.v  sinAh 

f  a"  ) 

Vi  V - 

7 r  . 

A  coshA/i 

1  <r  +  oj-  ' 

X 

(  (-  -  COSOlf 

-  —  sincof  )<(A . 

(23) 

to,"  ~  or  ==  f>k 

tanhA/i, 

(24) 

where  17  is  the  solution  found  by  Hammack  (1973) 
for  a  barotropic  ocean.  Hence,  the  surface  waves 
in  a  two-layer  system  are  not  affected  by  the  strati¬ 
fication  in  the  Boussinesq  limit.  Since  density  dif¬ 
ferences  are  typically  small  for  oceans  U  ~  0.002), 
previous  results  for  surface  waves  in  a  barotropic 
ocean  are  applicable  to  leading  order. 

It  is  also  useful  to  approximate  specific  properties 
of  the  internal  and  surface  waves  for  long  modes 
such  that 


e  a  =  O (A/t,)  =  O (A/ift  <§  1.  (25) 

P-.- 

Under  these  conditions  the  well-known  phase 
speeds  of  the  surface  and  internal  waves  are  found 
(from  Eq.  (22)]  to  be,  respectively. 


=  1 or  Ik-  -  git  1 


It 


—  +  CKe3) 
It, 


(26a) 


CV  =  ~ 


0(£-')|. 


(26b) 


The  significant  difference  in  phase  speeds  for  these 
long  waves.  CJC,  =  Ole1'4),  suggests  that  even 
though  the  surface  and  internal  waves  overlap 
spatially  near  the  generation  region,  they  will 
separate  1  ipidly  and  nonlinear  interactions  between 
(he  two  systems  will  not  have  time  to  become  sig¬ 
nificant 

In  addition  to  the  phase  speeds  it  is  also  of  interest 
to  examine  the  amplitude  of  the  motion  induced  by 
each  wave  system  at  the  opposite  interface.  Follow¬ 
ing  Keulegan  (1953),  normal-mode  analysis  sug¬ 
gests  that 

ft.,  =  —  .V, ,  (27a) 

It 


ft,  = 


€ 


where  0.,  is  the  coupled  motion  of  the  interface  due 
to  a  wave  /V,  at  the  free  surface  and  ft,  is  the  coupled 
motion  of  the  free  surface  due  to  an  internal  wave 
.Vj.  Here,  the  waves  ,V,  and  /V,  represent  the  inde¬ 
pendent  modes  at  the  free  surface  and  interface, 
respectively.  (Note  that  ,V,  =  17,  and  N2  =  r),  when 
the  surface  and  internal  waves  do  not  overlap 
spatially.)  The  motion  ft2  is  barotropic  and  would 
exist  at  the  interface  location  even  if  the  fluid  was 
not  stratified. 

In  order  to  see  the  wave  structure  in  ( 19),  numeri¬ 
cal  solutions  at  r  =  b  for  both  barotropic  and  baro- 
clinic  waves  propagating  from  the  source  region  are 
shown  in  Fig.  2.  Results  are  presented  for  impul¬ 
sive  response  at  the  free  surface  and  interface 
( Fig.  2a)  and  for  impulsive  interfacial  response  with 
creeping  free  surface  response  ( Fig.  2b).  Magnitudes 
of  the  nondimensional  generation  parameters  used  in 
computations  for  Fig.  2a  are  chosen  to  be  asymp¬ 
totically  representative  of  prototype  phenomena: 
dimensional  quantities  correspond  to  experiments 
which  will  be  discussed  in  Section  5. 

As  expected,  the  surface  waves  do  not  appear  to 
vary  with  layer-depth  ratio  in  Fig.  2.  In  fact  these 
results  are  identical  to  those  for  a  homogeneous 
ocean  computed  by  (23)  except  for  a  small  residual 
lowering  of  the  water  level  which  persis*s  until  the 
internal  wave  passes.  It  should  be  emphasized  that 
surface  and  internal  waves  overlap  spatially  in  the 
generation  region;  hence,  these  results  represent  a 
linear  combination  of  both  the  independent  (,V)  and 
induced  motions  (0)  at  the  interface  and  free 
surface.  The  residual  lowering  of  the  free  surface 
is  the  0(e)  motion  /3,  [see  Eq.  (27b)]  induced 
by  the  internal  wave.  In  fact,  this  induced  motion 
differs  for  each  depth  ratio;  however,  the  small 
differences  are  not  discernable  at  figure  scale.  The 
impact  of  the  surface  wave  at  the  interface, 
ft;,  is  more  pronounced  and  clearly  evident  in  Fig. 
2a.  This  induced  motion  is  easily  estimated  from 
(27a)  using  the  free  surface  solution  [since  iV,  =  77, 

+  O(e)].  When  ft2  is  eliminated  from  the  internal 
wave  computations,  iV2  results  and  is  also  shown  in 
Fig.  2a.  The  leading  waves  of  the  independent 
internal  and  surface  modes  are  similar  to  the  per¬ 
manent  deformation  of  the  sea  floor.  The  maximum 
amplitude  of  the  surface  wave  is  one-half  the  sea 
floor  uplift,  i.e.,  (,V,)inax  =  £„/2,  as  found  previously 
for  a  homogeneous  ocean.  The  maximum  amplitude 
of  the  independent  internal  waves  (/V .,)„,.„  for  these 
impulsive  motions  appears  to  vary  with  the  layer- 
depth  ratio  according  to 

GVftmax  =  <Mt)GVP> max  =  (£„/2)(  It  ,//l ).  (28) 


(27b) 
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This  important  result  and  further  clarification  of  the 
effects  produced  by  varying  the  size  scale  blh  as  well 
as  relative  thickness  of  the  upper  layer  may  be  estab¬ 
lished  in  the  following  manner.  Restricting  atten¬ 
tion  to  (mathematically)  impulsive  bed  motions, 
the  initial  ( t  =  0+)  amplitudes  at  x  =  b  can  be  eval¬ 
uated  in  closed  form  to  yield  in  the  Boussinesq 
limit 


2£„  ,  rr  b  i 

r)  i{h. O')  =  -  arctan  jtanh  —  —  J  , 


(29) 


T?,4/>.0*) 


Co 

=  —  arctan 


sinh(  77/>//t)/cos|  ~  ) 


(30) 


The  result  for  the  surface  wave  (29)  is  identical  to 
that  for  a  homogeneous  ocean  (Sells,  1965)  and 
clearly  shows  that  tj,  —  N,  —  17  —  Co/2  for  blh  >  1. 
Eqs.  (29)  and  (30)  are  shown  in  Fig.  3  as  a  function 
of  blh  and  hjh.  (Note  that  (29)  and  (30)  are  equiva¬ 
lent  when  h,lh  — >  0.]  Results  for  the  internal  wave 
are  identical  to  the  surface  wave  for  blh  >  1 
regardless  of  the  layer-depth  ratio.  For/;//;  <  1,  the 
internal  wave  amplitude  exceeds  that  at  the  free 
surface  as  hjh  increases  (e.g.,  hjh  >  0.25).  In  nil 
cases,  the  amplitudes  become  inversely  proportional 
to  the  size  scale  as  blh  decreases.  Noting  that 
iV -,(b. 0*)'  — T]2(b,0r)  -  f}2{b,0+),  where  p is  given 
by  (27a),  Eqs.  (29)  and  (30)  may  be  used  to  deter¬ 
mine  the  behavior  of  the  internal  mode  with  size 
scale;  the  result  is 


O') 

=  k  i 

•  n  j 

i 


arctan 


sinh(7r/;//i) 


cos 


| 

2  h  ) 


x  arctan 


(31) 


For  blh  >  1,  Eq.  (31)  yields,  in  closed  form,  the 
numerical  results  found  previously  for  the  internal 
wave  amplitude  given  by  (28).  Hence,  even  though 
the  free  surface  and  interface  are  lifted  (lowered) 
the  same  amount  by  an  impulsive  bed  motion  with 
blh  >  I.  part  of  the  potential  energy  gained  at  the 
interface  is  associated  with  the  barotropic  motion 
P<.  Consequently,  the  amplitude  of  the  internal 
wave  is  less  (by  the  factor  hjh)  than  that  of  the 
surface  wave  where  all  of  the  potential  energy 
gained,  to  O(e),  is  available  to  the  barotropic  mode. 


5.  Experimental  equipment  and  procedure 

A  seriesofexperiments  was  conducted  in  a  labora¬ 
tory  wave  tank  which  is  approximately  32  m  long, 
60  cm  deep  and  40  cm  wide  The  tank  is  equipped 
at  one  extremity  with  a  wave  generator  consisting  of 
a  movable  bed  section  (piston)  whose  time-dis¬ 
placement  history  is  controlled  by  an  electro- 
hydraulic-scrvo  system.  The  wave  generator  can 
accurately  model  the  bed  motion  described  by  (1). 
Both  the  wave  generator  and  tank  have  been 
described  in  detail  by  Hammack  (1972). 
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In  order  to  model  the  density  stratification  of  Fig. 
1,  freshwater  and  brine  are  used.  The  miscibility 
of  these  fluids  produces  a  finite  pycnocline  thick¬ 
ness  6;  thus,  this  stratification  provides  a  better  test 
of  the  two-layer  model  as  a  useful  approximation  for 
thin  pycnocline  regions.  The  tank  is  stratified  by 
first  introducing  a  layer  of  freshwater  with  a  depth 
/'1-  Brine  is  then  introduced  beneath  the  fresh¬ 
water  at  ports  distributed  along  the  tank  bed  until 
a  total  depth  h  is  attained.  The  volume  of  the 
reservoir  used  to  mix  the  brine  limits  the  depth  It, 
to  a  maximum  of  10  cm.  For  the  experiments 
presented  herein,  the  total  depth  h  is  fixed  at  15  cm 
while  the  depth  of  the  brine  is  varied  with  h<  =  5. 
7.5  and  10  cm.  A  density  difference  of  ~5%  is 
utilized  in  all  experiments.  This  difference  is  suf¬ 
ficiently  large  to  provide  stability  of  the  stratifica¬ 
tion  during  filling  yet  sufficientlysmall  to  provide  an 
adequate  test  of  analytical  results  based  on  small 
density  differences.  The  length  b  of  the  movable  bed 
section  is  constant  (61  cm)  for  all  experiments. 
Hence,  the  experimental  size  scale,  blh  =  4.07,  is 
sufficiently  large  for  results  to  be  typical  of  the 
tsunamigenic  case.  Density  stratifications  are  meas¬ 
ured  using  a  conductivity  probe  at  two  locations 
before  and  after  each  series  of  experiments;  typi¬ 
cal  results  are  presented  in  Fig.  4  where  /i,  =  5  cm 
and  li2  =  10  cm.  The  similarity  of  initial  stratifica¬ 
tions  at  both  locations  demonstrates  the  uniformity 
of  this  stratification  along  the  wave  tank.  A  thick¬ 
ness  8  for  the  pycnocline  region  is  calculated 
based  on  the  definition 


yielding  initial  values  of  80  =  0.9  and  1.0  cm  and 
final  values  of  8,  =2.1  and  1 .8  cm  for  the  upstream 
and  downstream  positions,  respectively.  Actual 
stratifications  are  qualitatively  similar  to  theoretical 
solutions  of  the  diffusion  equation  for  an  initial 
density  profile  with  a  simple  discontinuity  at  one 
level.  Using  this  solution  {p  =  p,  4-  dp  erfc[y/ 
(4Dt),  J J}  with  the  initial  discontinuity  occurring 
at  y  =  0,  the  characteristic  thickness  8  defined  by 
(32)  is  related  to  the  diffusion  coefficient  D  by 


8  =  (4-nDt)'  -.  (33) 

Hence,  an  apparent  diffusion  coefficient  can  be  cal¬ 
culated  for  the  time  interval  of  the  experiments  by 


I  /(S,)2  ~  (8„); 

4n  \  At 


(34) 
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The  apparent  diffusion  coefficients  for  the  measure¬ 
ments  of  Fig.  4  are  D  =  2.3  x  IO~->  cm2  s'*  at  the 
upstream  position  and  D  =  1.7  x  10-5  cm2  s1  at  the 
downstream  position.  The  more  intense  mixing  at 
x  =  b  results  from  the  shedding  of  a  vortex  at  the 
edge  of  the  bed  section  during  its  motion  (see 
Hammack,  1972).  The  mixing  of  the  interface  at  this 
position  was  visually  detected  by  adding  blue  dye  to 
the  brine.  At  the  downstream  position  the  apparent 
diffusion  coefficient  is  only  slightly  greater  than 
the  molecular  diffusion  coefficient  (D,„  =  1.5  x  10~5 
cm2  s'1)  for  salt  in  water.  Hence,  the  shear  layer 
developed  during  the  passage  of  internal  waves  is 
laminar  and  does  not  mix  the  fluids  appreciably. 
Each  series  of  experiments  was  stopped  when  the 
interfacial  thickness  at  x  =  b  reached  —2  cm. 

In  order  to  quantitatively  measure  the  internal 
waves,  a  laser-optics  detector  system  is  used.  A 
schematic  drawing  of  the  internal  wave  gage  is 
shown  in  Fig.  5.  The  light  beam  from  a  helium- 
neon  gas  laser  (wavelength  =  632.8  nm  with  0.5  mW 
output)  is  aligned  perpendicular  to  the  glass  side- 
walls  of  the  wave  tank.  A  system  of  cylindrical 
lenses  transforms  the  circular  beam  into  a  sheet  of 
light  with  a  constant  vertical  dimension  of  4.5  cm 
which  is  directed  horizontally  across  the  wave  tank. 
After  traversing  the  tank,  the  light  is  focused  onto 
a  photodiode  which  provides  an  output  voltage  pro¬ 
portional  to  the  incident  light  intensity.  By  dying  the 
brine  dark  blue  and  permitting  the  interface  to  inter¬ 
cept  the  light  sheet  when  the  fluids  are  quiescent, 
subsequent  motions  of  the  interface  induce  changes 
in  the  output  signal  of  the  photodiode  which  are 
displayed  on  an  oscillograph  recorder.  It  should  be 
noted  that  the  dyed  interface  seen  by  the  internal 
wave  gage  appears  in  the  upper  region  of  the  dif¬ 
fuse  salinity  interface  as  indicated  in  Fig.  4.  The 
laser,  optics  and  detector  system  are  mounted  on  a 
continuous  frame  shaped  in  an  inverted  U  over  the 
wave  tank;  hence,  the  wave  gage  is  calibrated  by 
raising  and  lowering  the  frame  differing  amounts 
during  quiescent  fluid  condition.  Calibration  curves 
are  inherently  nonlinear  for  this  measurement  sys¬ 
tem.  Internal  waves  with  maximum  heights  <  4.5 
cm  are  easily  measured;  however,  large-amplitude 
motions  are  generally  avoided  during  wave  measure¬ 
ments  at  t  =  b  due  to  the  vortex-induced  mixing 
at  the  position.  The  vortex  generates  a  high- 
frequency  signal  at  the  photodiode  which  is  super¬ 
posed  on  the  signal  corresponding  to  the  internal 
wave.  In  addition,  the  background  signal  for  quies¬ 
cent  conditions  may  change  appreciably  before  and 
after  an  experiment.  Surface  waves  are  electroni¬ 
cally  recorded  using  conventional  parallel-wire 
resistance  probes  and  an  oscillograph  recorder. 
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6.  Comparison  of  theoretical  and  experimental 
results 

In  order  to  examine  the  near-field  wave  struc¬ 
ture,  the  temporal  behavior  at  the  edge  of  the  gen¬ 
eration  region  (,r  =  b)  is  examined.  Wave  motions 
at  this  location  are  primarily  free  (unforced)  and 
provide  the  initial  condition  for  propagation  models 
which  may  be  patched  to  the  generation  model. 
Results  are  presented  for  a  wide  range  of  time-size 
ratios  in  order  to  more  clearly  delineate  the  tsunami- 
genic  case  (impulsive  with  It,  «  lu  and  bib  >  I). 

a.  Wave  profiles 

Comparisons  of  theoretical  wave  profiles  at 
v  =  b  with  experimental  data  are  shown  in  Figs.  6 
and  7.  Results  are  shown  for  positive  (uplift)  and 
negative  (downifyow)  bed  motions.  Linear  com¬ 
putations  for  the  free  surface  motion  agree  well  with 
measurements.  (Surface  wave  measurements  ate 
coincident  for  the  two  cases  of  bed  uplift.)  Linear 
results  for  the  internal  waves  also  agree  well  with 
measurements;  however,  there  are  differences 
which  require  explanation.  Measured  waves  at  the 
interface  exhibit  a  high-frequency  noise  of  small- 
amplitude  superposed  on  the  main  wave.  The  source 
of  this  high-frequency  noise  is  twofold.  First,  as 
mentioned  in  Section  5,  a  vortex  is  shed  at  the  down¬ 
stream  edge  of  the  piston  (.r  =  b)  during  rapid 
motions.  The  effect^of  this  vortex  on  the  amplitude 
measurement  at  the  interface  is  most  pronounced  for 
positive  motions  with  thin  lower  layers  (/i,  <  It,) 
and  was  visually  observed  to  be  responsible  for 
the  large  depression  in  amplitude  near Ug/h)' 2  =  12 
in  measurement  (d)  of  Fig.  6.  Second,  gas  bubbles 
formed  in  the  thin  gap  between  the  piston  and  its 
confining  tank  walls.  These  bubbles  were  freed 
during  bed  motion  and  subsequently  those  from  the 
gap  along  x  =  b  penetrated  the  light  sheet  of  the 
measurement  transducer.  (Brine  was  especially 
prone  to  bubble  formation  resulting  from  aeration 
during  mixing  and  chemical  reactions  at  unavoid 
able  scratches  in  the  anodized  aluminum  of  the 
piston.) 
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The  lime  required  for  the  lead  wave  to  propagate 
past  v  =  b  differs  from  the  linear  prediction  for 
both  surface  and  internal  waves  in  Fig.  6.  This 
behavior  was  observed  previously  for  surface  waves 
by  Hammack  (1973)  and  may  be  partially  attributed 
to  nonlinear  effects  arising  from  both  finite  wave 
amplitudes  and  finite  changes  in  depth  in  the  gen¬ 
eration  region.  Based  on  amplitude  considerations 
alone,  differences  between  linear  analysis  and 
measurements  are  expected  to  be  maximum  when 
/i,  S  h  ,  and  minimum  when  //,  =  h  ..  This  expected 
behavior  is  exhibited  in  Fig.  6  even  though  ampli¬ 
tude  nonlinearity  is  too  small  to  account  for  all  of  the 
observed  discrepancy.  Another  contributing  factor 
is  the  eddy  shed  at  x  =  h  during  piston  motion. 
Fluid  mixed  by  the  eddy  often  remains  in  this  re¬ 
gion  for  some  time  after  the  waves  had  departed. 
(Eventually,  this  fluid  spreads  into  its  layer  of  neutral 
buoyancy.)  It  was  confirmed  qualitatively  by  visual 
observation  that  this  phenomenon  causes  an  appar¬ 
ent  increase  in  both  the  wave  amplitude  near  its 
leeward  portion  and  the  lead  wave  period  similar  to 
the  behavior  observed  in  measurement  (</)  of  Fig. 
6  (where  /r ,  >  h2). 

Mode  separation  and  motions  induced  at  the  op¬ 
posite  interface  are  more  clearly  demonstrated  by 
the  downstream  (.v  =  b  +  20 h)  results  shown  in 
Fig.  8.  Both  experimental  measurements  and  linear 
computations  are  presented  for  an  impulsive  bed 
motion  with  £0  =  1.5  cm, /i,  =  5  cm  and  /i.,  =  10  cm. 
The  surface  wave  arrives  first  at  the  downstream 
station  at  t(glh)'  -  =  20  corresponding  to  the  speed 
C,  given  by  (26a);  the  induced  motion  f3,  at  the 
surface  is  clearly  evident.  Observed  and  predicted 
behavior  of  the  surface  wave  and  its  induced  mo¬ 
tion  agree  well.  Both  the  interface  and  free  surface 
become  quiescent  prior  to  the  arrival  of  the  in¬ 
ternal  wave.  Predicted  and  measured  shapes  of  the 
internal  wave  are  similar;  however,  the  maximum 
amplitude  of  the  measured  wave  is  20%  less  than 
predicted.  In  addition,  the  time  required  by  the 
observed  wave  to  reach  the  downstream  station 
corresponds  to  an  average  speed  which  is  10%  less 
than  the  linear  prediction  for  C2.  These  discrep¬ 
ancies  probably  result  from  the  viscous  and  non¬ 
linear  effects  present  in  the  experiments. 

h.  Maximum  amplitudes 

It  is  instructive  to  examine  the  variation  in  lead 
wave  amplitude  with  the  time-size  ratio  of  the  bed 
motion.  In  fact,  results  of  this  type  provide  the 
basis  for  classification  of  the  generation  process 
(impulsive,  creeping  and  transitional).  Fig.  9  shows 
experimental  and  theoretical  results  at  v  -  b  tor 
both  the  surface  and  internal  waves.  The  maximum 
amplitudes  of  the  independent  modes  (A/,  and  <V2) 
are  shown  as  a  function  of  their  respective  time- 
size  ratio  and  layer-depth  ratio.  Experimental  meas¬ 
urements  for  both  positive  and  negative  bed  motions 
are  presented.  (Nonlinear  effects  for  all  experiments 
shown  in  Fig.  9  are  expected  to  be  small  based 
on  the  discussion  in  Section  3.)  The  time-size 
ratio  for  the  surface  response  (r,)  spans  the 
classification  range  of  impulsive  to  creeping  while 
the  internal  response  (t2)  ranges  from  impulsive 

I  -'s'liu-'  il 
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Surface  wave  amplitudes  in  Fig.  9  are  identical  for 
each  layer-depth  ratio  (since  the  size  scale  bill  is 
constant)  and  equivalent  to  that  found  for  a  homo¬ 
geneous  ocean.  Both  experimental  and  analytical 
results  in  Fig.  9  demonstrate  that  CV,)nuis  =  £„/ 2  for 
impulsive  motions  while  ( .V,  )max  a.  £0/r,  in  the 
creeping  regime;  previous  results  indicate  that  this 
behavior  i?  independent  of  size  scale  when  bib  >  1 . 
The  theoretical  and  experimental  results  for  the 
internal  wave  amplitudes  agree  well.  In  addition,  the 
variation  of  (/Vj)mitx/£n  with  r-<  is  similar  to  that  of  the 
surface  waves  except  for  an  amplitude  reduction  by 
the  factor  hjh  which  appears  to  be  valid  even  for 
time-size  ratios  exceeding  the  impulsive  regime. 
This  behavior  permits  the  barotropic  and  baroclinic 
results  to  be  collapsed  into  a  single  curve  as  shown 
in  Fig.  10  for  the  full  range  of  r,  and  r2  (when  the 
size  scale  satisfies  bill  >  1).  The  weak  dependence 
(narrowing)  of  the  transitional  regime  with  larger 
size  scales  C blh  1)  is  also  indicated  in  Fig.  10. 

7.  A  comment  on  Coriolis  effects  and  the  two-layer 
approximation 

It  is  customary  and  appropriate  to  neglect  Coriolis 
effects  for  barotropic  tsunamis  since  their  energy  is 
concentrated  at  periods  (say.  /,  =  30-60  min)  small 
compared  to  possible  inertial  periods.  (Alterna¬ 
tively,  the  length  scale  L  of  the  barotropic  mode 
is  small  compared  to  the  Kelvin-Rossby  radius  of 
deformation.)  However,  the  sluggishness  of  the 
baroclinic  response  indicated  by  C2/C,  —  e1'*  =  0.04 
immediately  suggests  the  possibility  of  mesoscale 
periods  for  these  modes.  Since  the  length  scale  L 
for  both  barotropic  and  baroclinic  modes  is  the 
same  and  fixed  by  the  source  mechanism,  a  first 
estimate  of  the  dominant  baroclinic  period  is 

h  ~~  ——  1 ,  ~  «-"*/,  '  25/,. 
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When /,  is  of  order  or  greater  than  one  inertial  period 
U,  =  12  h  (sin<l>r'  where  <t>  is  the  source  latitude] 
the  baroclinic  modes  must  be  considered  as  inertio- 
gravity  waves.  In  fact,  Coriolis  effects  will  limit  the 
period  of  these  inertio-gravity  waves  to  t ...  <  t, 
with  t.,  — »  i,  as  L  — »  (e.g.,  see  Kraus,  1972,  p. 

215).  Neither  the  results  of  the  preceding  analysis 
nor  the  waveguide  (one-dimensional  propagation) 
approximation  are  applicable  for  describing  the 
long-term  (an  inertial  period  on  greater)  dynamics 
ot  these  waves.  However,  the  preceding  results 
still  remains  applicable  for  propagation  times  much 
less  than  an  inertial  period.  In  either  case,  the 
variaodity  of  both  source  lengths  for  tsunamigenic 
earthquakes  and  the  inertial  period  (with  latitude) 
suggest  that  the  potential  role  of  Coriolis  effects 
is  earthquake  and  site  specific. 

Finally,  it  should  be  noted  that  the  two-layer 
stratification  adopted  herein  permits  an  unbounded 
frequency  spectrum  for  the  baroclinic  modes.  In 
actuality,  the  continuous  stratification  of  oceans 
will  limit  free  baroclinic  oscillations  to  frequencies 
less  than  the  maximum  Brunt-Vaisala  frequency, 
say  N,„.  [As  noted  by  Kraus  (1972,  p.  218),  the  two- 
layer  model  is  equivalent  to  the  first  baroclinic 
mode  of  a  continuous  stratification  with  equivalent 
phase  speeds.]  For  a  typical  ocean  stratification 
with  ,V,„  =  0.01  s  ',  baroclinic  wave  energy  will  be 
restricted  to  periods  c,  >  10  min.  As  already  im¬ 
plied,  energy  generally  will  be  concentrated  at  pe¬ 
riods  much  greater  than  this  lower  bound. 
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Fig.  I.  Fluid  domain  and  sea  floor  dislocation  model: 

(a)  spatial,  (b)  temporal. 

Flu.  2.  Theoretical  wave  profiles  at  x  =  b  =  61  cm  for «  =  0.05  and  h  =  15  cm: 

( - 1  free  surface  waves  (all  It, lit):  broken  curves  represent  internal  waves: 

l - )  =  2/3:  ( - )  It, lit  =  1/2:  ( - )  Mi  =  1/3:  I . )  V..  (a) 

t ,  =  0.1.  t.  =  0.01:  (b)  r,  =  5.2,  t,  =  0.5. 

Fig.  3.  Theoretical  Variation  of  maximum  wave  amplitudes  with  source  size 
scale  hilt  at  v  =  h  for  impulsive  response:  (solid  line),  (tj,  )ma,/{o'.  dashed  line. 
( 1 1 ■: 

Fig.  4.  Experimental  density  profiles  for  h,  ~  5  cm  and  h.  =  10  cm:  line  with 
open  circles,  initial  profile:  line  with  closed  circles,  final  profile  (a)  v  =  h.  St 
=  205  mjp;  (b)  t  =  b  +  120/i,  St  -  185  min.  S  indicates  elevation  of  ob¬ 
served  (dyed)  interface 


Fig.  5.  Schematic  drawing  of  internal  wave  gage 

Fig  6.  Theoretical  and  experimental  waves  at  x  =  6  =  61  cm 
with  h  =  15  cm  and  e  =  0.05:  solid  line,  theory:  dashed  line, 
experiment,  (a)  Free  surface;  (b)  interface:  /i ,  =  5  cm.  C„  =  1.00 
cm.  r,  =  0. 13.  r,  -  0.0(4.  (c)  interface:  h,  =  7.5  cm.  Cn  =  -  I  00 
cm.  r,  =  0.25.  t,  =  0.027;  (d)  interface:  h,  -  10  cm.  Co  =  1.08 
cm,  r,  =  0.12.  r,  =  0.012. 

Fig  7.  Theoretical  and  experimental  waves  at  i  =  h  =  61  cm 
with  h  -  15  cm  and  e  =  0.05:  solid  line,  theory,  dashed  line, 
experiment,  (a)  Free  surface;  (b)  interface  It,  -  5  cm.  Cn  =  -LOO 
cm,  t,  =  5.13.  r2  =  0.53,  (c)  interface,  h,  =  7.5  cm,  Co  =  0.97 

cm.  r,  =  4.53,  Tj  =  0.50;  (d)  interface:  It,  =  10  cm,  Co  =  2.00 
cm.  r,  =  4  87,  r,  =  0.50. 
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Fig.  K.  Theoretical  and  experimental  waves  at  x  -  b  +  20;i: 

It  =  15  cm.  ft,  =  5  cm.  h  =  61  cm,  c  —  0.05,  U  =  1.5  cm. 
r,  =  0  12.  t,  =  0.012;  solid  line,  free  surface:  dashed  line,  inter¬ 
face  tat  Linear  theory:  (b)  experiment. 

Fit.  9  Variation  of  with  r,  and  with  t_.  at  i  =  b  =  61  cm  for  It  -  15  cm. 

t  -  0  05.  Solid  curve  represents  surface  wave  calculation  (for  all  /ip/n  using  linear  theory  Vertical 
slash  indicates  experimental  data  for  surface  waves 

FtO.  10.  Theoretical  variation  of  ( .V , ).na, /£„  with  r,  and  l(<V})milv/(„)(/iiKi,|  with  t;.  at  t  =  h.  Dashed 
curve  indicates  asvmntotic  results  for  h  li  t>  1 
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1.  Introduction 


Studies  of  tsunami  generation  often  employ  simple  models  of  the  sea 
floor  dislocations  to  permit  tractable  analytical  solutions.  Although 
these  solutions  provide  basic  insight  into  the  generation  process,  they 
are  incapable  of  producing  explicit  results  for  prototype  events  where 
both  the  spatial  and  temporal  distributions  of  the  sea  floor  dislocation 
may  be  quite  complicated.  Herein  we  exploit  the  apparent  linearity  of 
the  generation  process  and  demonstrate  both  the  use  and  validity  of  the 
superposition  principle  to  construct  solutions  for  complex  bed  motions. 
Analytical  and  experimental  results  are  presented  for  a  monopolar  dis¬ 
location  (block  upthrust  or  downthrow)  with  a  complex  time-displacement 
history.  The  time  history  used  in  the  computations  is  obtained  from  an 
integrated  accelerogram  recorded  at  Pacoima  Dam,  near  Los  Angeles, 
during  the  earthquake  of  February  9,  1971.  A  complex  spatial  deforma¬ 
tion  is  not  used  in  order  to  enable  experimental  verification  of  the 
analytical  results.  This  is  unfortunate  since  it  appears  that  the 
details  of  the  time-displacement  history  are  not  important  for  proto¬ 
type  phenomena  where  the  motion  may  be  considered  instantaneous.  How¬ 
ever,  it  is  important  to  note  that  the  analysis  treats  both  space  and 
time  variations  in  an  identical  manner;  hence,  confirmation  of  this 
approach  for  complex  time  variations  strongly  suggests  analogous  behav¬ 
ior  for  complex  spatial  variations. 

Finally,  we  examine  and  compare  several  alternative  time-displacement 
histories  for  the  mean  motion.  It  is  shown  that  the  results  for  each  mean 
motion  can  be  unified  by  introducing  a  velocity  as  a  descriptive  param¬ 
eter  which  is  based  on  the  kinetic  energy  input  of  the  moving  bottom  to 
the  overlying  fluid. 


2.  Time-Displacement  History 

To  define  a  complicated  time-displacement  history  with  relevance  to 
earthquake- induced  ground  motions,  we  have  chosen  an  accelerogram  for 
the  vertical  component  of  motion  recorded  at  Pacoima  Dam  (near  Los 
Angeles)  during  the  earthquake  of  February  9,  1971.  The  accelerogram 
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component  of  ground  acceleration  recorded  2/9/71  0600  PST  Pacoima  Dam 
a,  with  integrated  velocity  and  ground  displacement. 


shown  in  Figure  1  may  be  integrated  numerically  using  appropriate  pre¬ 
cautions  ( e.g .  ,  see  Nigam  and  Jennings,  1968),  to  yield  a  time  history 
of  the  vertical  ground  displacement.  Results  of  the  integration  pro¬ 
cess  for  the  "smoothed"  velocity  history  and  the  corresponding  ground 
displacement  history  are  shown  below  the  accelerogram  in  Figure  1. 

An  analytical  representation  of  the  algebraically  complex  time- 
displacement  history  of  Figure  1  over  a  finite  record  interval, 

0  <  t  £  T,  may  be  constructed  using  a  Fourier  series  of  the  form: 

N 

eft)  =  C  +  l  C  sin  (K  t  +  3  ).  (1) 

o  ,  n  n  n 

n=l 

In  (1)  Cn  are  the  amplitudes  of  the  Fourier  components  with  CQ  repre¬ 
senting  the  mean  (permanent)  ground  displacment,  Kn  are  the  component 
wave  frequencies,  and  3^  are  the  component  phase  angles.  The  accuracy 
of  (1)  in  representing  the  integrated  displacement  of  Figure  1  is 
determined  by  the  number  of  components  N  retained  in  the  Fourier  sum. 
Results  of  computations  with  N  =  18  are  illustrated  in  Figure  2;  this 
truncated  sum  will  be  adopted  in  the  subsequent  analysis. 

It  should  be  emphasized  that  both  the  instrument  characteristics 
and  the  numerical  integration  techniques  used  to  obtain  ground  dis¬ 
placements  from  acceleration  measurements  necessarily  distort  (filter) 
information  in  long  period  components.  In  particular,  the  mean  (and 
permanent)  displacement  of  the  integrated  motion  shown  in  Figure  1  (and, 
of  course,  its  Fourier  representation  in  Figure  2)  is  not  expected  to 
accurately  model  the  actual  permanent  deformation.  To  compensate  for 
this  distortion,  we  may  again  exploit  the  superposition  principle  and 
add  a  nonsinusoidal  component  to  the  Fourier  series  representation  of 
(1).  As  an  example  of  this  approach,  consider  the  ramp  motion  in  time 
of  the  form: 


Cr(t)  =  C0t/T  (2) 

during  the  time  interval  0  <  t  £  T;  other  choices  for  the  mean  motion 
are  discussed  and  compared  in  section  5.  In  practice,  the  sum  of  the 
mean  component  resulting  from  the  Fourier  synthesis,  C0,  and  the  addi¬ 
tional  component  Q0  from  (2)  should  be  chosen  to  equal  the  actual  per¬ 
manent  ground  offset.  Hence,  a  general  representation  for  a  complicated 
time-displacement  history  becomes: 


C(t) 


c  + 

o 


C0t/T 


•  £. 

n=l 


C  sin 
n 


(K  t 
n 


(3) 


3.  Solutions  of  the  Water  Wave  Problem 


Consider  a  two-dimensional  (x,y)  and  incompressible  ocean  of  uniform 
depth  h  initially  in  equilibrium  with  the  earth's  gravitational  field  g 
which  acts  in  the  negative  y  direction.  At  time  t  =  0  a  section  of  the 


isplacement  ( - )  and  Fourie 


sea  floor  begins  to  deform  vertically  with  a  time  and  spatial  variation 
given  by  £(x,t).  We  seek  the  inviscid,  irrotational ,  and  barotropic 
deviations  ri(x,t)  of  the  ocean  free  surface  from  its  equilibrium  posi¬ 
tion.  With  the  coordinate  system  at  the  equilibrium  position  of  the 
free  surface,  the  linearized  description  of  motion  in  terms  of  a  veloc¬ 
ity  potential  <p  =  4>(x,y,t)  is: 


4>xx(x,y,t)  +  4>yy(x,y,t)  =  o 

(4) 

4>y(x,-h,t)  =  Ct(x,t) 

(5) 

4>y(x,0,t)  =  nt(x,t) 

(6) 

4>t(x,0,t)  =  -gn(x,t) 

(?) 

where  subscripted  variables  indicate  partial  differentiation, 
convenient  to  eliminate  n(x,t)  in  (6)  and  (7)  by  combining  to 
single  condition  for  the  velocity  potential: 

It  is 
yield  a 

4>tt(x,0,t)  +  g<J)y  (x,0,t)  =  0 

(8) 

Using  the  Laplace  transform  in  t  and  the  Fourier  transform  in 
tions  (4),  (5)  and  (8)  become: 

x,  equa- 

?yy(k,y,s)  -  k^(k,y,s)  =  0 

(9) 

^y  (k  ,-h ,  s)  =  sc;  (k  ,  s) 

(10) 

^y(k,0,s)  +  <jT(k,0,s)  =  0 

(11) 

where  the  overbar  of  a  function  f(x,t)  indicates: 

f(k,s)  =  |  dx  |°°  e^X  e  St  f(x,t)dt. 

(12) 

-"O  -oo 


Solving  (9),  (10),  and  (11)  for  <J>(k,y,sl  and  noting  from  (7)  that 

n(k,s)  =  -  (s/g)^(k,0,s)  (13) 

we  find: 

t 

r)(k,s)  =  s2  C(k,s)/(s2  +  u)2)  cosh  kh  (14) 


where  w2  =  gk  tanh  kh.  Inverting  the  Laplace  and  Fourier  transforms 
yields: 


n(x,t) 


r  co  r  . 

_L  f  /lim 
2 it  J 

_oo 


1 

2fri 


|U+ir 

M-ir 


s2  e~lkx  est  qk,s)  ds| 
(s2  +  a)2)  cosh  kh  i 


dk 


(15) 
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Again,  the  final  integration  over  wavenumber  k  must  be  evaluated  numeri¬ 
cally. 


4,  Comparison  of  Theory  and  Experiment 

As  noted  earlier,  the  simple  block  deformation  of  the  sea  floor  was 
chosen  to  enable  experimental  verification  of  the  analytical  model 
developed  for  complicated  time-displacement  histories.  The  experimental 
facility  used  in  these  tests  has  been  described  in  detail  by  Hammack 
(1972,  1973).  Basically,  the  wave-maker  consists  of  a  rectangular 
piston  in  the  bottom  of  a  wave  tank  (and  spanning  its  width)  whose 
motion  is  controlled  by  an  electro-hydraulic-servo  system.  The  servo 
system  converts  a  time-voltage  command  signal  into  a  proportional  verti¬ 
cal  displacement  of  the  piston.  For  the  experiments  reported  herein, 
the  piston  length  in  the  direction  of  wave  motion  is  b  =  61  cm  while  the 
quiescent  water  depth  above  the  piston  is  h  =  10  cm.  Before  presenting 
results  of  the  tests,  we  describe  the  motivation  for  choosing  other 
experimental  scales. 

A  "global"  time  scale  for  the  forcing  of  the  overlying  ocean  by  the 
sea  floor  is  the  period  T.  The  appropriate  time  scale  for  the  gravita¬ 
tional  response  of  the  long  barotropic  wave  modes  is  b/fgh)^  which  cor¬ 
responds  to  the  time  required  for  waves  to  escape  the  generation 
region.  For  prototypical  earthquakes  the  ratio  of  the  forcing  and 
response  time  scales,  t  =  Tfgh^/b,  termed  the  time-size  ratio,  is  small 
so  that  details  of  the  time-displacement  history  generally  are  not 
important.  However,  our  interests  herein  require  that  the  details  of 
the  temporal  motion  have  a  significant  impact  on  the  generated  wave 
structure.  Hence,  the  period  T  for  the  experimental  tests  must  be 
scaled  so  that  t  exceeds  unity;  in  fact,  for  the  experiments  a  period 
T  =  4  secs  was  chosen  which  yields  t  =  6.5.  Previous  experiments  by 
Hammack  (1973)  also  indicate  that  the  generation  process  for  proto¬ 
typical  tsunamis  is  linear  and  that  nonlinearity  remains  insignificant 
for  vertical  displacements  which  do  not  exceed  about  20%  of  the  over- 
lying  ocean  depth.  This  criterion  is  adhered  to  in  the  experiments  by 
restricting  the  instantaneous  displacement  of  the  piston  to  less  than 
2  cm. 

In  the  first  test  we  examine  experimental  and  theoretical  results 
for  the  time-displacement  history  shown  in  Figure  2  using  experimental 
parameters  T  =  4  secs  and  a  chosen  permanent  (mean)  displacement  of 
CQ  =  0.38  cm,  with  the  amplitudes  and  frequencies  Kj,  of  the  eighteen 
Fourier  components  used  in  Figure  2  scaled  appropriately.  (In  Figure  2 
the  corresponding  parameters  of  the  actual  ground  displacement  are 
T  =  40  secs  and  CQ  =  7.63  cm.)  The  scaled  Fourier  components  are  then 
summed  and  the  result  is  converted  to  an  analogue  (time-voltage) 
signal  which  is  used  to  command  the  wavemaker.  The  resulting  wave 
motion  at  the  leading  edge  of  the  piston  (x  =  b)  is  measured,  and  the 
results  are  shown  in  Figure  3.  Theoretical  results  at  x  =  b  are 
evaluated  for  each  of  the  eighteen  (scaled)  Fourier  components  accord¬ 
ing  to  (17)  and  summed  to  yield 


the  results  also  are  shown  in  Figure  3.  The  excellent  agreement  between 
the  predicted  and  measured  data  is  self-evident. 

Similar  results  are  shown  in  Figure  4  where  a  ramp  mean  motion  with 
an  amplitude  of  Co  =  1.33  cm  is  added  to  the  Fourier  synthesis  of  Figure 
2.  The  total  permanent  displacement  of  the  piston  is  C o’  =  +  Co  = 

1.71  cm  which  has  been  used  to  normalize  the  measured  and  theoretical 
wave  amplitudes.  The  theoretical  result  is  equivalent  to  (25)  with  an 
added  component  for  the  ramp  computed  from  (24).  Again,  the  agreement 
between  measured  and  computed  data  is  excellent  with  the  wave  structure 
clearly  showing  the  added  volume  (mass)  resulting  from  the  enhanced  mean 
displacement. 


5.  A  Comment  on  Mean  Motions 

In  previous  studies  (Hammack,  1972,  1973)  two  additional  models  for 
the  mean  displacement  of  a  block  section  of  the  sea  floor  have  been 
examined.  These  time  histories  are: 

a.  exponential:  Ce(t)  =  CQ[1  -  exp  ( - 1 . 1 t/Te) ] 

b.  half-sine:  Cg(t)  =  CQ[(1  -  cos  irt/Ts)H(Ts-t)/2  +  H(t-Ts)] 

and  we  repeat  for  completeness  the  mean  motion  introduced  here: 

c.  ramp:  C/t)  =  CQ [tH(Tr-t)/Tr  -  H(t-Tr)]. 

The  three  mean  motions  listed  above  span  a  wide  range  of  displacement 
characteristics.  We  note  that  the  choice  of  characteristic  time  scales 
Te,  Ts  and  T  are,  in  fact,  arbitrary  to  a  certain  extent  even  though 
"natural"  choices  are  apparent.  (This  flexibility  is  most  obvious  for 
the  exponential  motion  where  Te  was  chosen  for  experimental  convenience 
to  represent  the  time  for  two-thirds  of  the  displacement  to  occur.)  It 
has  been  found  that  the  properties  of  waves  generated  by  these  motions 
correlated  strongly  with  the  time-size  ratio  x  based  on  these  time 
scales.  For  example,  the  maximum  wave  amplitude,  say  q0,  occurring  at 
x  =  b  when  normalized  by  the  permanent  displacement  £0  exhibits  a  simple 
and  similar  functional  dependence  on  x  for  each  bed  motion.  With  the 
size  scale  b/h  of  the  dislocation  fixed  and  for  x<<1,  the  normalized 
amplitude,  r\0/C,0,  teaches  a  maximum  value  of  one-half  for  all  size  scales 
exceeding  unity.  Bed  motions  with  x<<1  are  termed  impulsive.  For  x 
very  large,  termed  creeping  generation,  Ho/^o  decreases  at  a  rate  which 
is  inversely  proportional  to  x.  The  constant  of  proportionality  for 
creeping  generation  varies  with  the  specific  choice  of  the  characteris¬ 
tic  time  scale  for  the  mean  motion.  Since  it  is  unlikely  that  any  of 
these  mean  motion  models  is  "correct"  from  a  geophysical  point  of  view, 
there  is  a  need  to  seek  a  unification  of  results  by  generalizing  the 
concept  of  characteristic  time  scale.  One  generalization  which  closely 


where  the  complex  inversion  integral  for  the  Laplace  transform  has  been 
used.  In  (15)  ui  is  the  wave  frequency  (w  =  is)  and  k  is  the  wavenumber. 
Explicit  results  for  specific  deformations  of  the  sea  floor  £(x,t)  will 
be  developed  now. 

3.1  Solution  for  a  single  Fourier  component 

Consider  a  block  section  of  the  sea  floor  of  length  2b  whose  time- 
displacement  history  corresponds  to  a  single  Fourier  component  of  (1). 
With  the  coordinate  system  centered  above  the  block  section,  we  have 

Cn(x,t)  =  CnH(b2-  x2)[sin(KT+6n)H(T- t)  +  sinf^t  +  0n)H(t  -  T)  ]  (16) 

where  H(  )  is  the  Heaviside  step  function.  Finding  the  transform  of 
(16) ,  substituting  into  (15)  ,  performing  the  integration  around  the 
Bromwich  contour,  taking  only  the  real  part  of  the  resulting  integral, 
and  noting  that  the  integrand  is  an  even  function  of  k,  we  find: 

%<*■«  ■  ^  f  °°fllrh1farb  {**B-H(t-T)[C.D.E]}  dk  (17) 

o 

where 
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sin  Bn [ (to2  sin  ait  - 

K  2  sin 
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Knt)/(W2  -  Kn2)] 

(18) 

B  = 

K  cosg  [(a)  sin  out  - 
n  n 

kn  sin 
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(19) 

C  = 

sln(KT  .  6n)j 

[[u2  cos  to(t-T)  -  Kn2  cos  Kn(t-T)]/(U)2-Kn2  ) 

}  (20) 

D  = 

Kncos(lt„™n>i 
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[gj  sin 

oj(t-T) 

-  K  sin  K  (t-T)  ]  /  (uj2-  K  2  )  1 
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(21) 

E  = 

-  sin  (K  T  +  8  )  cos 
n  n' 

u(t-T) 

. 

(22) 

The  final  integration  over  wavenumber  k  in  (17)  is  obtained  by  numerical 
quadratures. 

3.2  Solution  for  ramp  mean  motion 

The  ramp  time-displacement  history  of  (2)  for  the  block  deformation 
is  described  by 

C(x,t)  =  H(b2  -  x2)  [tH(T  -  t)/T  ♦  H (t  -  T)].  (23) 

Following  the  same  procedure  outline  in  section  3.1,  we  obtain  (again 
after  considerable  algebra)  the  water  surface  motion  nr  due  to  the  ramp 

V-o  ■  ~ir  Prlos-ra-  (Jrji51"  “<t-T»dk-'24) 

o 
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RESULTANT  WATER  WAVE  AT  x  = 


Figure  3.  Theoretical  ( - )  and  experimental  ( — )  wave  records  at 


Theoretical  ( - )  and  experimental  (-— )  wave  records  at  x  =  b  with  ramp  mean  motion. 


Theoretical  and  experimental  variati 


produces  the  desired  unification  is  the  following.  Consider  a  velocity 
V,  where  V2  is  the  kinetic  energy  imparted  to  the  fluid  by  the  sea  floor 
divided  by  one-half  the  total  mass  of  fluid  displaced  during  bed  motion 
PspbC0);  then  V  is  given  by: 

1 

v  =  ~  J  Ct2  dC  (26) 

and  a  corresponding  time  scale  T  =  C0/V.  In  terms  of  this  time  scale 
the  time-size  ratio  becomes 

T*  =  T(gh)%/b  =  C0(gh)i5/bV  (27) 

The  variation  of  p0/C0  at  x  =  b  for  b/h  =  6.1  with  t*  is  shown  in  Figure 
S  for  each  of  the  mean  motions;  both  theoretical  and  experimental 
results  are  presented.  For  all  of  the  experimental  data  we  have  taken 
£0/h  <  0.2  to  avoid  significant  nonlinear  effects.  The  collapse  of 
results  for  such  a  wide  range  of  mean  motion  characteristics  shown  in 
Figure  5  is  good  although  a  small  spread  still  exists. 


6.  Conclusions 


We  have  demonstrated  both  the  application  and  validity  of  a  strategy 
which  employs  multiple  uses  of  the  superposition  principle  to  develop 
theoretical  solutions  for  waves  generated  by  sea  floor  motions  with  com¬ 
plicated  time-displacement  histories.  Although  a  more  useful  test  for 
prototypical  phenomena  would  utilize  complicated  spatial  distributions 
for  the  sea  floor  dislocation,  the  tests  herein  were  restricted  to 
simple  block  dislocations  due  to  experimental  limitations.  However,  it 
is  emphasized  that  the  solution  method  does  not  distinguish  between 
space  and  time,  and  the  validity  established  herein  strongly  suggests 
that  the  methods  could  be  extended  to  complicated  spatial  deformations. 
Finally,  we  have  demonstrated  that  wave  properties  (in  particular  the 
maximum  amplitude  of  waves  escaping  the  generation  region)  for  a  wide 
range  of  mean  motion  characteristics  may  be  (almost)  collapsed  into  a 
single  functional  relationship  in  terms  of  a  time-size  ratio  based  on  an 
average  vertical  velocity  of  the  sea  floor  obtained  from  energy  consid¬ 
erations. 
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International  School  of  Physics 
"Enrico  Fermi" 

Course  LXXX,  Summer  1980 
Topics  in  Ocean  Physics 

SMALL-SCALE  OCEAN  WAVES 

J.  L.  Hammack 

Department  of  Civil  Engineering,  University  of  California  -  Berkeley,  California 

I.  —  Introduction. 

Presented  below  are  three  lectures  on  various  aspects  and  types  of  ocean 
waves.  In  all  cases  we  are  concerned  with  small-scale  waves  for  which  gravi¬ 
tation  provides  the  dominant  restoring  force.  The  relevant  length  and  time 
scales  of  these  waves  are  small  compared  to  the  Kelvin-Fossby  radius  of  defor¬ 
mation  and  inertial  period,  respectively,  and  Coriolis  effects  may  be 
neglected.  In  addition,  the  small  length  scales  (relative  to  Earth's  radius) 
permit  us  to  ignore  the  curvature  of  the  ocean  surface  and  adopt  a  plane- 
earth  approximation.  Various  other  approximations  will  be  adopted  as 
appropriate  in  order  to  obtain  model  equations  which  are  tractable  analytically. 
(Fortunately,  these  tractable  models  appear  to  remain  relevant  for  geophysical 
phenomena!)  In  particular,  we  are  interested  in  weakly  nonlinear  systems  in 
which  the  nonlinear  effects  manifest  slowly;  both  nonlinear  self-interacting 
and  resonant- triad  systems  are  considered. 

An  outline  of  these  lectures  is  as  follows.  In  Part  II  we  examine  two 
nonlinear  model  equations  for  the  evolution  of  gravity  waves:  the  Korteweg- 
deVries  (KdV)  equation  for  long  barotropic  and  baroclinic  waves  and  the 
nonlinear  Schrodinger  (NLS)  equation  for  short  barotropic  waves.  After  a 
brief  review  of  the  asymptotic  (large  time)  solution  of  these  equations  by 
inverse  scattering  theory,  experimental  data  are  presented  which  demonstrate 
the  reality  of  soliton  predictions — at  least  on  laboratory  scales.  In  Part  III 


2 


we  exploit  the  exact  solution  of  the  KdV  equation  and  its  linear  approximations 
in  order  to  derive  rather  precise  criteria  for  modelling  the  evolution  of 
long-wave  initial  data  on  geophysical  or  laboratory  scales.  The  modelling 
criteria  of  Part  II  are  then  applied  to  a  typical  oceanic  tsunami  in  order 
to  choose  relevant  model  equations  for  propagating  the  (barotropic)  wave  from 
its  generation  region  to  specific  target  sites.  Finally,  in  Part  IV  we 
examine  the  nonlinear  excitation  of  "edge"  waves  nearshore  by  linear  wave 
trains  incident  from  deep  water.  Both  theoretical  and  experimental  results 
are  presented  which  document  the  excitation  of  two  progressive  edge-wave  inodes 
through  a  nonlinear  resonant-triad  interaction  with  the  wave  reflected  after 
normal  incidence  from  offshore. 

Much  of  the  material  in  Parts  II  and  III  has  appeared  previously  in 
the  literature.  Hence,  many  of  the  details  concerning  experimental  equipment, 
procedures,  and  analysis  is  omitted;  a  thorough  discussion  may  be  found  in 
the  cited  references.  The  material  presented  in  Part  IV  on  edge  waves  is 
recent  and  not  presently  available  in  the  literature.  For  clarity,  a  more 
detailed  discussion  of  this  material  is  presented.  We  also  note  here  that 
the  notation  between  Parts  II,  III,  and  IV  is  not  necessarily  consistent. 

II.  — Water  Wave  Solitons. 

1.  —  Long  Gravity  Waves 

Consider  two  fluid  layers  with  uniform  densities  resting  in  a 

gravitationally  stable  configuration  (see  Figure  1)  on  a  horizontal  and 
impermeable  bed  of  infinite  lateral  extent.  The  upper  layer  possesses  a 
free  surface  along  which  pressure  is  constant;  surface  energy  effects  on 


3  - 


and  the  .fluid- fluid  interface  are  negligible.  Required  are  the  two- 

dimensional,  inviscid,  irrotational  motions  which  are  bounded  everywhere  and 

evolve  from  given  initial  distribution  of  velocities,  free  surface  deformation, 

and  interfacial  deformation — the  classical  water  wave  problem.  In  particular, 

we  are  interested  in  the  barotropic  displacement  of  from  its  static 

equilibrium  position  and  the  baroclinic  displacement  r\  of  from  its 

equilibrium  position.  To  derive  the  KdV  equation  as  the  appropriate  model 

for  both  displacements  from  the  governing  equations,  the  following  assumptions 

are  necessary.  First,  characteristic  wave  lengths  k  ^  (an  inverse  wave 

2 

number)  must  be  long  relative  to  the  total  fluid  depth;  i.e.,  (kh)  <<  1, 

so  that  dispersive  effects  are  weak.  Second,  characteristic  wave  amplitudes 

a  are  small  relative  to  the  total  depth;  i.e.,  a/h  <<  1,  so  that  nonlinearity 

is  weak.  Third,  both  weak  effects  of  dispersion  and  nonlinearity  are  approxi- 

2 

mately  equal;  i.e.,  e  =  a/h  'v  (kh)  .  For  simplicity,  we  will  further  assume 

at  the  outset  that  density  differences  in  our  two-layer  ocean  are  small  (as  in 

its  geophysical  counterpart)  so  that  A  =  (p^  -  <<  1-  Expanding  the 

dependent  parameters  in  the  governing  equations  in  a  power  series  in  e  and 

introducing  the  multiple  time  scales  tg  *  t,  t  =»  et,  .  .  . ,  one  finds  at 

leading  order  that  the  evolution  equations  for  the  initial  wave  are  hyperbolic 

(nondispersive)  and  linear.  At  this  order,  which  corresponds  to  the  fast  time 

scale  tg,  an  initial  disturbance  decomposes  into  four  modes  consisting  of 

left-  and  right-running,  barotropic  and  baroclinic  modes.  The  phase  speeds 

C,  and  C  of  these  modes  are 
1  2 


2 

(la)  Barotropic  (surface)  waves:  =  gh 

2 

(lb)  Baroclinic  (interfacial)  waves:  =  g'h^h^/h 


where  we  have  invoked  the  Boussinesq  limit  A  -*■  0  with  gA  =  g*  <<  g  remaining 

finite.  All  wave  modes  propagate  with  permanent  form  and  do  not  interact 

with  each  other  or  themselves;  the  baroclinic  inodes  propagate  much  slower, 

1/2 

0(e  ),  than  the  barotropic  modes. 

2 

At  the  next  order  (s  )  weak  nonlinear  effects  and  dispersion  occur. 

Each  wave  mode  experiences  a  self-interaction  on  the  slow  time  scale  t  =  et, 
but  no  interactions  between  modes  occur  due  to  their  rapid  separation  by  the 
phase  speed  differences  of  (1) .  The  self-interaction  of  the  right-running 
barotropic  mode  is  governed  by  a  dimensional  equation  of  the  form 


(2) 


ni  +  cini  +  +  hh\ 

t  X  X  XXX 


=  0 


A  more  convenient  choice  of  nondimens ional  variables  for  describing  these 


waves  are 


<3>  < 


X  =  (x  -  C  t)/h,  r  =  7(g/h)1'/2t 

J-  6 


f(x»T)  =  jn(x,t)/h 


V 

In  terms  of  these  variables,  (2)  reduces  to  the  KdV  equation  with  the  common 
form 


(4)  f  +  6ff  +  f  =  0 

t  X  XXX 


In  a  similar  manner,  the  dimensional  equation  governing  the  self- 

i 


interaction  of  the  right-running  baroclinic  waves  is 


n_  +  c  n  +  -c  (- - —  )n_n~  +  tCJi  h  n 

2^  22  22  h_  h,  22  62122 

t  x  2  1  x  xxx 


Introducing  the  normalized  variables 


x  =  (x  -  c2t)/(h1h2)‘,‘ 

1  1/2 

(6)  <  t  =  -Hg'/hJ^t 
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,  3  1  1  . 

f  =  -(- - 7— )n 

2  h2  hj 

equation  (5)  reduces  to  the  KdV  equation  given  by  (4) . 

The  asymptotic  solution  of  the  KdV  equation  for  arbitrary  initial  data 
f(X/0)  H  fQ  by  inverse  scattering  theory  has  been  described  in  detail  by 
Segur'*’  inter  alios.  Here,  we  briefly  list  features  of  the  asymptotic  solution 
to  be  illustrated  in  the  laboratory  experiments. 

i)  An  initial  disturbance  evolves  into  a  finite  number  of  permanent 

waves  (solitons)  ordered  by  their  amplitude.  When  the  solitons  are 
well  separated,  the  local  shape  of  each  is  given  by 

2  1/2 

(7)  f  =  asech  {(a/2)  (X  “  XQ  -  2ax) } 

where  a  and  Xq  are  constants.  The  rank  ordered  solitons  are 
followed  by  a  dispersive  train  (radiation)  of  oscillatory  waves, 
ii)  The  number  N  of  solitons  evolving  from  initial  data  of  finite 

extent,  say  f  =  0  for  x  <  X^  and  x  >  X2'  i-s  equivalent  to  the  number 
of  zeros  of  $  for  X  >  X^  where  <J>  satisfies 


±A  + 
dx2 


2  +  VXH  =  ° 


l(xj»  lM  (xJ  =  0. 

L  dv  1 


iii)  When  the  net  volume  V  (or  mass)  in  the  initial  wave  is  finite  and 


positive,  i.e. 


(9) 


V  = 


fQ(x)dx  >  0 


at  least  one  soliton  emerges. 

iv)  When  f  <  0  everywhere,  no  solitons  emerge  and  the  asymptotic 
solution  consists  (only)  of  the  radiation  components, 

v)  Two  other  important  classes  of  data  are  those  for  which  V  <  0  with 
f q 4^  0  for  all  x  and  those  for  which  V  =  0.  No  general  statements 
regarding  the  asymptotic  solution  for  these  cases  is  provided 
theoretically.  Experiments  demonstrate  that  the  evolution  of 
solitons  depends  on  the  detailed  structure  of  the  initial  data. 


1.1  Experimental  Procedures  In  order  to  illustrate  the  applicability 

cf  the  KdV  equation  as  a  model  for  long  water  waves,  a  series  of  experiments 

is  conducted  in  a  laboratory  wave  tank  31.6  m  long,  61  cm  deep,  and  39.4  cm 

wide.  For  studying  the  barotropic  motions  in  a  system  with  small  density 

differences,  the  density  stratification  plays  no  role;  hence,  in  these 

experiments  a  uniform  density  fluid  of  depth  h  is  used.  A  detailed  descrip- 

2 

tion  of  these  experiments  is  given  by  Hammack  and  Segur.  In  order  to  study 
the  evolution  of  baroclinic  waves,  the  tank  is  stratified  with  fresh  water 
and  brine.  Details  of  these  experiments  can  be  found  in  Hammack^  as  well  as 

4 

Hammack  and  Segur. 

In  the  long  wave  experiments  reported  here,  the  wavemaker  consists  of 
a  rectangular  piston  located  in  the  tank  bottom  adjacent  to  an  upstream 
end-wall.  The  piston  spans  the  tank  width  and  has  a  length  b  in  the  direction 


of  wave  motion  (see  Figure  2) ;  lengths  of  b  =  30.5  cm  and  61  cm  are  used. 

The  vertical  motion  of  the  piston  is  controlled  by  am  electro-hydraulic-servo 
system  and  completely  user  specified.  In  a  typical  experiment,  the  piston 
is  moved  for  a  finite  time  interval.  Differences  between  the  initial  and 
final  position  of  the  piston  permit  the  net  volume  V  in  the  generated  wave 
train  to  be  calculated.  Various  initial  data  are  generated  by  varying  the 
time-displacement  history  of  the  piston. 

Both  surface  and  internal  wave  amplitudes  are  measured  at  fixed  loca¬ 
tions  (stations)  along  the  tank.  We  note  here  that  differences  between  the 
temporal  variation  in  wave  amplitude  at  a  fixed  spatial  position  and  the 
spatial  variation  and  wave  amplitude  at  a  corresponding  fixed  time  are  small, 
0(e),  and  neglected  in  all  calculations. 

1.2  Results  for  Barotropic  Wave  Evolution 

1.2.1  Initial  data  with  V  >  0. 

Figures  3  and  4  illustrate  the  evolution  of  two  barotropic  waves  with 
a  net  positive  volume.  Normalized  wave  amplitudes  are  presented  at  four 
succeeding  stations  along  the  tank  in  a  coordinate  system  which  moves  with 
the  linear  (nondispersive)  speed  C^.  Note  that  the  leading  portion  of  the 
wave  system  appears  at  the  left  in  these  figures.  Shifts  of  the  waves  to  the 
left  (right)  at  succeeding  stations  indicate  phase  speeds  greater  (less)  than 
in  this  coordinate  system. 

The  initial  wave  at  (x  -  b)/h  =  0  in  Figure  3  is  rectangular  and  positive 
and  appears  to  sort  itself  into  three  positive  waves  identified  by  separate 
crests  (local  maxima)  during  propagation.  These  three  waves  are  rank  ordered 
by  amplitude  and  are  followed  by  a  weak  train  of  dispersive  waves  whose  speed 
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is  much  less  than  C^.  In  order  to  examine  the  local  shape  of  the  leading  two 
waves,  theoretical  soliton  profiles  defined  by  (7)  have  been  superposed  on 
the  measured  data  at  the  last  measurement  station  (using  the  measured  wave 
amplitudes  for  specifying  the  parameter  a) .  These  profiles  are  not  extended 
into  regions  where  the  two  waves  are  still  interacting  strongly  with  adjacent 
wave  structure.  Clearly,  the  measured  waves  appear  to  be  locally  KdV  soli tons. 
Further  evidence  that  the  three  lead  waves  at  the  last  measurement  station 
are  solitons  is  provided  by  numerically  integrating  (8)  using  the  wave  profile 
at  (x  -  b)/h  =  0  as  the  initial  data  f^.  Computations  indicate  that  N  =  3 
solitons  should  evolve  in  agreement  with  the  observed  pattern  of  evolution. 

The  evolution  of  a  more  complicated  initial  wave  with  V  >  0  is  shown 
in  Figure  4  where  a  leading  negative  wave  is  followed  by  a  (larger)  positive 
wave.  After  only  twenty  depths  of  propagation,  the  positive  wave  has 
separated  into  three  separate  crests  while  the  negative  wave  remains  essen¬ 
tially  unchanged.  During  subsequent  propagation  the  three  labelled  crests 
of  the  positive  wave  appear  to  retain  their  integrity  as  they  progress  thro’ 
the  leading  negative  wave  and  emerge  at  the  front  of  the  wave  train.  At  the 
last  measurement  station,  labelled  waves  1  and  2  clearly  resemble  KdV 
solitons.  The  third  wave  is  still  interacting  with  the  once-leading  negative 
wave  at  the  last  station  and  cannot  be  unmistakably  identified  as  a  soliton. 
However,  computations  with  (8)  using  the  wave  at  (x  -  b)/h  =  0  as  f^  yield 
N  =  3  strongly  suggesting  the  third  wave  is  indeed  a  soliton.  We  note  that 
other  experiments  also  indicate  that  solitons  evolve  from  the  positive  waves 
in  the  initial  wave  and  can  be  identified  long  before  asymptotic  conditions 


fmmm  - 


are  achieved. 


1.2.2  Initial  waves  with  V  <  0 


Results  for  the  evolution  of  two  barotropic  wave  systems  with  V  <  0  are 
shown  in  Figures  5  and  6.  The  initial  wave  in  Figure  5  is  the  negative  counter¬ 
part  of  the  experiment  shown  in  Figure  3  where  three  solitons  appeared  to 
evolve.  Over  similar  distances  of  propagation,  no  solitons  appear  to  evolve 
in  Figure  5  —  just  as  expected  for  initial  data  with  f  <  0  for  all  x- 
Instead,  a  negative  wave  evolves  whose  frontal  slope  decreases  and  lengthens 
with  time.  This  lead  wave  is  followed  by  a  train  of  strongly  dispersive 
waves  with  phase  speeds  much  less  than  C^.  In  fact,  the  wave  structure  of 
Figure  5  represents  the  radiation  solution  of  the  KdV  equation  as  shown  with 
quantitative  tests  by  Hammack  and  Segur. 5 

Further  evidence  for  the  evolution  of  complicated  waves  with  V  <  0  is 
shown  in  Figure  6.  In  this  case  the  lead  positive  wave  evolves  into  one  or 
possibly  two  waves  while  the  negative  wave  evolves  in  a  manner  similar  to 
that  observed  in  Figure  5.  Computations  using  the  initial  wave  in  (8)  indicate 
that  one  soliton  should  develop  asymptotically.  Hence,  the  second  labeled 
wave  is  expected  to  eventually  disappear  as  it  appears  to  be  doing. 

1.2.3  Initial  waves  with  V  =  0 

Figure  7  illustrates  the  evolution  of  initial  data  with  V  =  0  consisting 
of  a  positive  wave  preceded  and  followed  by  negative  waves.  Applications  of 
(8)  to  the  initial  wave  suggest  that  two  solitons  should  evolve.  The  large 
positive  wave  in  the  initial  data  quickly  separates  into  two  crests  which 
appear  to  migrate  through  the  lead  negative  wave  during  subsequent  propagation. 
At  the  last  station,  labeled  wave  1  has  progressed  to  the  front  and  clearly 
has  the  shape  of  a  soliton  everywhere.  The  second  labeled  wave  still  appears 


10 


to  be  interacting  with  the  frontal  slope  of  the  once  leading  negative  wave. 

A  trailing  grain  of  dispersive  waves  similar  to  that  of  Figure  5  also  evolves. 

1.3  Results  for  Baroclinic  Wave  Evolution  Ihe  evolution  of  baroclinic 

long  waves  at  the  interface  of  a  stratified  fluid  with  h^  =  45  cm,  h^  =  5  cm, 

and  A  =  0.05  is  shown  in  Figure  8.  Wave  amplitudes  recorded  at  seven 

stations  along  the  tank  in  a  coordinate  system  that  moves  with  the  linear 

speed  are  presented.  The  initial  wave  is  positive,  and  according  to 

calculations  by  (8),  should  evolve  two  solitons.  Ihe  observed  evolution  is 

in  agreement  with  this  prediction.  (A  strongly  damped  train  of  oscillatory 

waves  has  been  omitted  in  Figure  8.)  A  more  quantitative  comparison  of  the 

lead  wave  profile  with  the  theoretical  shape  (7)  at  the  last  four  stations  of 

measurements  is  shown  in  Figure  9.  The  agreement  with  the  measured  data  in 

Figure  9  is  excellent.  Further  results  for  baroclinic  wave  evolution  are 

4 

presented  in  Hammack  and  Segur. 

2.  —  Short  Gravity  Waves 

We  now  turn  our  attention  to  the  evolution  of  short  gravity  waves;  i.e., 
waves  whose  characteristic  length  k  ^  is  comparable  to  or  greater  than  the 
local  ocean  depth;  only  barotropic  modes  are  considered.  In  fact,  our 
interest  in  solitons  requires  kh  >  1.36  which  we  adopt  as  the  definition  of 
the  short-wave  regime  herein.  In  order  to  derive  the  nonlinear  Schrodinger 
(NLS)  equation  as  an  evolution  model,  we  consider  a  wave  train  with  a 
dominant  and  identifiable  (mean)  wave  number  k^.  As  for  long  waves,  we  require 
weak  nonlinearity,  which  is  now  characterized  by  the  wave  steepness  e  =  akQ  <<  1. 
The  wave  system  is  permitted  to  have  weak  modulations  such  that  the  variation 
in  wave  number  6k  is  small;  i.e.,  5kA0  <<:  1*  In  other  words,,  we  are  concerned 


with  "narrow-band"  wave  systems  that  are  weakly  nonlinear.  Thirdly,  we 

postulate  a  balance  of  both  small  effects  so  that  e  =  ak^  *\»  5k/kg.  Multiple 

scale  analysis  of  the  governing  equations  again  yields  a  heirarchy  of 

problems  at  different  orders  of  e.  At  lowest  order  we  recover  the  linear 

£ 

dispersive  waves  of  Stokes  with  an  amplitude  a  that  is  constant.  At  the 

next  time  scale,  t  =  st,  the  wave  amplitude  is  modulated,  and  we  find  that 

the  amplitude  modulations  propagate  with  the  linear  group  speed  C  . 

2 

Continuing  to  the  third  order  with  a  time  scale  t^  =  e  t,  one  finds  that  the 
complex  amplitude  modulation  a  must  satisfy  the  NLS  equation.  Defining  non- 
dimensional  coordinates  as 

f 

X  =  ek0(x  -  Cgt)  . 

(10)  /  t  =  e2(gk0)1/2t 

A  =  k02(gk0r1/2a 

V 

the  NLS  equation  takes  the  form 

(11)  iA  +  AA  +v|a|  A  =  0 

XX 

where  A  and  v  are  known  functions  of  the  water  depth  h,  gravitation  g,  and 
carrier  wave  number  kQ  (or  frequency  u)Q).  Details  of  the  derivation  for 
finite  depth  can  be  found  in  Hasimoto  and  Ono^  inter  alios. 

Like  the  KdV  equation,  the  NLS  equation  cam  be  solved  exactly  for 
arbitrary  initial  data  by  inverse  scattering  theory.  When  kQh  >  1.36, 


envelope  solitons  can  occur;  the  one-soliton  solution  of  (11)  is 
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2\  '  2 
(12)  A  =  a  —  sech{ax}exp{iAa  t} 

where  a  is  an  arbitrary  constant  related  to  the  envelope  amplitude.  In  general, 
initial  data  of  finite  extent  will  evolve  a  finite  number  of  envelope  solitons, 
rank  ordered  by  the  group  velocities  of  their  dominant  carrier  waves  and 
embedded  in  a  dispersive  train  (radiation)  of  oscillatory  waves  which  decays 
in  amplitude  with  time.  It  is  important  to  note  that  the  speeds  of  both  the 
solitons  and  the  radiation  components  are  not  impacted  by  nonlinearity  at  this 
order.  Hence,  unlike  their  long  wave  counterparts  of  the  KdV  equation,  these 
solitons  and  radiation  components  do  not  separate  with  time.  This  suggests 
that  in  order  to  observe  clearly  the  evolution  of  envelope  solitons  even  in 
a  contrived  laboratory  experiment,  one  must  design  the  initial  data  such  that 
the  wave  content  of  the  radiation  spectrum  at  the  dominant  frequency  of  the 
soliton  carrier  waves  is  small.  Alternatively,  one  must  observe  evolution 
until  the  inviscid  decay  of  the  radiation  (by  frequency  dispersion)  combined 
with  viscous  decay  has  progressed  sufficiently. 

Although  a  quantitative  comparison  of  theory  and  experiment  for  the 

NLS  equation  analygous  to  that  of  the  KdV  equation  has  not  been  performed, 

0 

qualitative  tests  by  Yuen  et  al  provide  evidence  for  its  applicability  to 
narrow  band  systems.  We  present  in  Figure  10  an  (unpublished)  experiment 
which  illustrates  the  evolution  of  an  initial  wave  packet  whose  envelope 
amplitude  and  carrier  wave  frequency  are  modulated;  the  mean  wave  frequency 
is  =  0.6  sec.  The  experiment  is  conducted  in  a  tank  approximately 
50  m  long,  2.5  m  wide,  and  a  water  depth  of  h  =  1  m.  After  45  m  of  propaga¬ 
tion,  the  initial  packet  evolves  into  a  collection  of  wave  groups  which  are 
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ordered  by  the  group  speed  of  their  dominant  carrier  wave.  Energy  is  concen¬ 
trated  in  the  group  with  carrier  waves  at  the  dominant  period  u>  1  *  0.6  sec  of 
the  initial  data.  The  envelope  of  this  group  clearly  has  the  shape  of  the 
soliton  profile  (12)  shown  superposed  on  the  measured  wave.  The  envelopes 
for  the  other  wave  groups  do  not  agree  with  (12)  ,  and  apparently  these  groups 
represent  the  radiation  components  in  the  initial  data. 

Although  testing  of  the  NLS  equation  is  not  complete,  there  is  growing 

9 

evidence  (see  Yuen  and  Lake  )  that  the  dynamics  of  short-water  waves  is  much 
more  complicated  than  that  of  long  waves, and  more  complicated  models  are 
required.  The  restriction  of  the  NLS  equation  to  narrow  band  systems  and  its 
probable)  instability  to  transverse  perturbations  appear  especially  severe  for 
many  geophysical  applications. 

XXI.  —  Modelling  Criteria  for  Long  Water  Waves. 

Even  with  all  the  assumptions  implicit  in  the  classical  water-wave 
problem  introduced  in  Part  II,  the  general  equations  remain  intractable  analyt¬ 
ically.  Hence,  we  resort  to  further  approximations  such  as  those  required  to 
yield  the  KdV  and  NLS  equations.  Generally,  approximations  are  formalized  by 
perturbation  expansions  in  terms  of  a  small  parameter (s) .  In  application  of 
these  approximate  models,  questions  naturally  arise  as  to  when  the  inequalities 
used  in  ordering  the  physics  are  actually  satisfied;  i.e.,  how  small  is  small? 
For  evolution  models,  the  relevant  question  distills  into  "during  what  time 
interval  does  a  particular  approximation  correctly  represent  the  general 
solution?"  Closely  related  to  this  topic  are  questions  of  similitude  between 


laboratory  models  and  their  geophysical -scale  counterparts.  The  answers  to  all 
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of  these  questions  involve  determination  of  nondimensional  parameters  which 
characterize  the  phenomenon  in  question.  In  practice  these  parameters  are 
often  deduced  through  the  ingenious  use  of  dimensional  analysis,  by  examina¬ 
tion  of  the  governing  equations,  or  using  the  most  preferred  but  least  available 
method,  examination  of  exact  solutions  of  the  general  and  approximate  model 
equations.  Here  we  examine  the  KdV  equation  and  exploit  exact  solutions  of  its 
various  approximations  in  order  to  develop  rather  precise  criteria  for 
modelling  the  propagation  of  long  water  waves.  These  criteria  are  then  applied 
to  the  problem  of  tsunami  propagation  across  ocean  basins — the  topic  which 


precipitated  the  analysis.  Further  details  are  presented  by  Hammack  and  Segur.' 
1.  —  Analysis. 


In  order  to  obtain  definitive  results,  we  adopt  the  following  point  of 

view.  If  the  initial  wave  f  is  sufficiently  smooth  and  localized,  and  if 

2 

=  a/h  and  =  (kh)  based  on  the  initial  data  are  both  small,  then  the 
KdV  equation  (4)  will  be  the  approximate  model  eventually.  (In  fact,  we  know 
that  (4)  is  appropriate  when  e  =  0(e ^).)  If  and  e are  not  the  same 

order  of  magnitude  initially,  then  simpler  forms  of  the  KdV  equation  may  be 
applicable  for  some  time  interval,  e.g.. 


f  =  0  (linear,  nondispersive  model) 


f  +  f 
I  XXX 

f  +  ff 
I  X 


0  (linear,  dispersive  model) 

0  (nonlinear,  nondispersive  model) 


Ursell  showed  that  the  nondimensional  parameter 


U  =  =  ak  2/h3 
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which  we  shall  call  the  Ursell  number,  is  an  important  indicator  to  determine 
which  of  these  models  is  appropriate.  One  should  use  (14)  ,  (4) ,  or  (15) , 
depending  on  whether  U  <<  1,  U  =  0(1),  or  U  >>  1,  respectively.  Ursell  also 
provided  evidence  that  U  is  time  dependent  and  will  tend  to  an  order-one 
limit  so  that  all  waves  in  this  category  eventually  propagate  according  to 
(4)  as  postulated  previously.  Although  there  is  general  agreement  with 
Ursell' s  results,  there  is  disagreement  on  how  to  interpret  "order  unity"  and 
how  to  define  the  relevant  length  scale  k  ^  for  the  evolving  waves.  To  make 
his  criteria  more  precise,  we  examine  here  initial  data  for  which  <<  <<  1 

initially  so  that  the  linearized  models  (13)  and  (14)  are  the  relevant  approx¬ 
imations  of  (4) .  The  dimensional  length  scale  kQ  1  for  the  initial  data  is 
defined  as  its  overall  length,  and  the  dimensional  amplitude  scale  is  aQ. 
Further,  we  will  focus  on  criteria  for  modelling  the  leading  wave  only. 

Results  for  the  trailing  wave  structure  and  the  nonlinear  approximation  (15) 
are  given  in  Hammack  and  Segur.10 

1.1  Linear  Dispersive  Theory  Since  nonlinearity  is  small  for  the 
postulated  initial  data,  we  assume  a  small  parameter  y  «  1  exists  and  seek  a 
formal  series  solution  of  the  KdV  equation  (4)  in  the  form 

(17)  f  =  yfL  +  u2f2  +  0 (y3) 

(The  required  definition  of  y  will  come  out  of  the  analysis.)  Substituting 
(17)  into  (4)  yields  a  heirarchy  of  problems. 


(18) 

(f,  )  + 

(f,) 

1  X 

1  XXX 

(19) 

(f,)  + 

(f,) 

2  T 

2  XXX 

etc. 

=  0,yf1(x,0)  =  fQ 

=  -6W  ,f2(X.O)  = 


0 
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The  solution  of  (18),  which  is  equivalent  to  (14),  is  well  known: 

CO 

(20)  UfI  =  27  {  f0  (K)  exP  [*■  (<X  +  <3T)]dK 

—CO 

whe  re 

OO 

(21)  fQ(<)  =  {  fQ(x)exp(-iKx)dx 

—  OO 

is  the  Fourier  transform  of  the  initial  data.  For  practical  reasons,  we  are 
most  interested  in  the  asymptotic  form  (t  -*•  “>)  of  the  solutions  for  both  the 
linear  and  nonlinear  models.  As  t  •>  »  with  (x|/t  ->-0  in  order  to  remain  at  the 
wave  front,  the  asymptotic  form  of  (20)  is 


(22)  MfL(x»T)  =  fQ(0) (3T)‘1/3Ai(C)  -  if^(0) (3T)"2/3Ai' (C) 

-  |f^(0)  (3T)_1Ai"(S)  +  0  [  (3t)  ~4/3] 

1/3 

where  C  =  X/(3t)  and  Ai(£)is  the  Airy  function.  The  coefficients  in  (22) 
have  simple  interpretations: 


(23a) 


(23b) 


(23c) 


where  V  is  the  nondimensional  volume  of  the  initial  wave,  is  an  Ursell 
number  for  the  initial  wave,  and  8^,  8^  are  constants  that  depend  on  the 
details  of  the  initial  wave.  Thus  (22)  becomes 


(24)  Uf1(X.T)  =  V(3T)~1/3{Ai(C)  +  61(^r)  (3T)"1/3Ai'  (?) 

+  (3T)“2/3Ai"(C)  +  .  .  .} 

From  (24)  it  follows  that  the  time  required  for  this  representation  to  become 
asymptotic  (second  term  smaller  than  first,  etc.)  is  at  least 


(3t)1/3  »  UQ/|v| 


A  particular  solution  of  (19)  is 


!f2  (x»f)  =  -  {  /  pf^U/Tjdz}^ 


and  defining 


1  ,x 

Q(x)  «  ~  { /  f0u)dZr 


V  —  oo 


the  homogeneous  solution  of  (19)  is 

X  2  »  „ 

(28)  U2f  (X.t)  =  V2  /  Ai  (z)  dz  +  J~f  (0(k)  ~  (i*c>_1> 


exp  [i  (icy  +  <  t)  ]  dK 


Thus,  as  t  +  »,  jx|/T  0 


U2f  (X,t)  =  V2{-  [  /  Ai (z) dz] 2  +  /  Ai(z)dz  + 


Comparing  (24)  and  (29)  ,  the  appropriate  definition  of  y  is  found  to  be  simply 


V  =  V 


In  other  words,  linear  dispersive  theory  along  with  its  large-time  asymptotics 
requires  the  dimensionless  wave  volume  V  to  be  small.  (Note  that  V  is  independent 
of  time.)  Since  the  solution  of  (29)  remains  0(1)  as  t  -*■  <*>,  the  series  (17) 
cannot  remain  asymptotic  after 

(31)  (3t)1/3  %  | v| -1 

Thus,  if  the  initial  wave  has  a  dimensionless  volume  V  (assumed  ^  0)  and 
Ursell  number  Uq  based  on  its  initial  (given)  dimensions,  asymptotic  linear 
dispersive  theory  is  valid  in  an  interval  no  longer  than 

(32)  Uq/M  «  (3t)  1/3  «  l/|v|. 

1.2  Nondispersive  Linear  Theory.  Since  our  initial  wave  is  postulated 
to  be  linear,  then  UQ  <<  1  and  necessarily  [ V f  «<  1,  as  well.  The  derivation 
of  the  KdV  equation  outlined  in  Part  II  indicated  that  linear  nondispersive 
theory  (13)  occurs  on  the  first  time  scale  of  evolution.  Hence  we  may  write 

(33)  3t  «  e2  =  (kQh)  2  =  (UQ/V)2 

for  (13)  to  be  applicable. 

2.  Summary  and  Application  of  Criteria  to  Tsunamic  Propagation 

In  summary,  we  have  examined  the  evolution  of  long-wave  initial  data  which 
is  parameterized  by  an  initial  Ursell  number  UQ  <<  1  and  volume  (or  mass) 

^  <<<:  may  model  the  evolution  of  this  initial  data  using  linear  nondispersive 

theory  (f  =0)  during  a  time  interval 


-  l9  - 


0  <  3t  <<  Ti 


The  next  relevant  model  is  linear  dispersive  theory  (f  +  f  =01 

t  XXX  ’ 


which  is  valid  during 


T i  <<  3x  <<  T2  =  U0  3 


with  its  asymptotics  becoming  valid  during 


t2  <<  3t  <<  t ^  =  y2. 


Subsequently,  we  must  use  the  KdV  equation  (f  +  6ff  +  f  =0)  for 

t  X  XXX 


3x  »  T  - 


(The  question  of  the  time  scale  for  applicability  of  KdV  asymptotics  is  discussed 

by  Hammack  and  Segur^.)  As  an  example  of  the  application  of  these  results, 

.  .  .  12 
consider  the  dimensional  scales  adopted  by  Carrier  for  major  tsunamis  which 

4-1  5 

impact  entire  ocean  basins:  aQ  =  10  ft. ,  h  =  1. 5  '  10  ft.,  k  =  2  *  10  ft. 

Then  UQ  t  0.1  and  V  %  0.01  and  linear  nondispersive  theory  is  valid  for  times 

2 

corresponding  to  propagation  distances  (using  =  gh)  of 


0  <  x  <<  600  miles. 


Linear  dispersive  theory  is  appropriate  for 


600  miles  <<  x  <<  6,000  miles 


indicating  that  dispersion  may  affect  this  wave  over  much  of  typical  ocean 
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trajectories.  Even  so,  linear  asymptotics  do  not  apply  until  x  >>  6,000  miles, 
which  exceeds  the  length  of  realistic  trajectories.  Hence,  linear  dispersive 
asymptotics  and  the  KdV  equation  are  not  required  for  describing  the  lead  wave 
of  the  tsunami  discussed  here.  (The  KdV  equation  may  be  used,  but  it  is 
unnecessarily  complicated.)  Of  course,  other  factors  may  be  required  to 
accurately  model  long-term  tsunami  propagation  (such  as  the  variable  bathymetry 
along  trajectories  and  three-dimensional  spreading  of  wave  energy) ;  the  intent 
of  the  analysis  here  is  only  to  develop  insight  into  the  relative  importance 
of  dispersion  and  nonlinearity. 

IV.  —  Excitation  of  Standing  Edge  Waves  on  Beaches 

1.  —  Introduction 

In  recent  years  considerable  attention  has  been  focused  on  the  occurrence 
of  ocean  waves  which  become  trapped  and  capable  of  concentrating  energy  in 
localized  regions  (wave  guides) .  Wave  trapping  can  occur  whenever  gradients 
exist  in  a  parameter  which  affects  the  wave's  phase  speed,  e.g.,  Coriolis  parameter, 
Brunt- Vaisala  frequency,  current  speeds,  and  water  depth.  Stokes  provided 
the  first  theoretical  evidence  of  trapping  for  surface  (gravity)  waves  near  the 
shoreline  of  a  plane  sloping  beach.  Stokes  found  a  normal  mode  solution  for  the 
(barotropic)  departure  n  of  the  water  surface  from  its  static  equilibrium  position 
(see  figure  11a)  of  the  form 

(38)  n  (x,y,t)  =  a  sin  6  exp  (-k  y  cos  8)  cos  k  x  sin  w  t 

m  m  mm 

where  -tan  8  is  the  beach  slope,  x  points  in  the  longshore  direction,  and 
y  points  offshore.  The  (linear)  dispersion  relation  for  the  various  longshore 
mode  number  m  =  1,2,  ...  is 


-  19 


(34) 


0  <  3t  <<  t i 


fn  1 


The  next  relevant  model  is  linear  dispersive  theory  (f  +  f  =  0) 

t  XXX  ' 

which  is  valid  during 


(35) 


T  1  <<  3t  <<  7  2  = 


V 


with  its  asymptotics  becoming  valid  during 


(36) 


T2  <<  3x  <<  T. 


V' 


Subsequently,  we  must  use  the  KdV  equation  (f  +  6ff  +  f  =0)  for 
"  t  x  XXX 


(37) 


3x  >>  T. 


(The  question  of  the  time  scale  for  applicability  of  KdV  asymptotics  is  discussed 
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by  Hammack  and  Segur  .)  As  an  example  of  the  application  of  these  results, 
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consider  the  dimensional  scales  adopted  by  Carrier  for  major  tsunamis  which 


4-1  5 

impact  entire  ocean  basins:  aQ  =  10  ft.,  h  =  1.5  '  10  ft.,  kQ  =2-10  ft. 


Then  Uq  'v  0.1  and  V  ^  0.01  and  linear  nondispersive  theory  is  valid  for  times 


corresponding  to  propagation  distances  (using  =  gh)  of 


0  <  x  <<  600  miles. 


Linear  dispersive  theory  is  appropriate  for 
600  miles  <<  x  <<  6,000  miles 

indicating  that  dispersion  may  affect  this  wave  over  much  of  typical  ocean 


2 

(39)  to  =  gk  sin  B 

m  m 

where  k  =  imr/b  is  fixed  by  the  beach  width  b.  Note  that  the  Stokes'  mode 
m 

of  (38)  is  periodic  in  the  longshore  direction  with  crests  pointing  offshore 

which  decay  in  amplitude  with  an  e- folding  distance  y  =  (k  cos  B)--*-. 

e  m 

The  crest  amplitude  is  maximum  at  the  shoreline  with  a  magnitude  a  sin  B 

where  a  is  referred  to  as  the  run-up  amplitude.  (Of  course,  run-up  phenomena 

cannot  be  represented  by  linear  solutions;  however,  the  reality  of  run-up  is 

acknowledged  since  a  is  the  up-beach  length  of  the  horizontal  projection  of  the 
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vertical  shoreline  amplitude.)  Ursell  demonstrated  that  the  Stokes’  mode 
is  only  the  lowest  (j  =  0)  of  a  discrete  set  j  =  -0,1  ...  ,  J  of  trapped  modes 
where  J  is  the  greatest  integer  satisfying 

IT 

(40)  (2J  +  1)  6  <  2  . 

These  higher  discrete  modes  are  algebraically  complicated  but  retain  longshore 

periodic  behavior.  Their  crest  amplitude  is  maximum  at  the  shoreline  but 

oscillates  in  the  offshore  direction  with  j  nodes  while  decaying  exponentially 

(see  figure  lib).  To  complete  the  set  of  normal  modes,  Hanson”'-'’  showed  that 

there  exists  a  continuous  spectrum  of  waves  with  wave-number  magnitude  k 
2 

satisfying  w  >  gk.  These  modes  are  even  more  complicated  to  describe 
algebraically  but  resemble  simple  deep-water  wave  trains  far  offshore.  An 
example  of  one  of  these  modes  with  crests  parallel  to  the  shoreline  is  shown 
in  figure  lib.  We  emphasize  here  that  the  continuous  spectrum  modes  are  not 
trapped  like  the  discre'te  modes . 
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Oceanographic  interest  in  edge  waves  appears  to  have  originated  in 

Isaacs  Zt  aZ. 16  suggestion  that  these  waves  might  be  responsible  for  "surf  beat" 

phenomena.  Since  then,  edge  waves  have  been  indicted  as  a  potential  mechanism 

for  numerous  features  of  coastal  dynamics  such  as  beach  cusps  (Guza  and  Inman)17, 

crescentic  bars  (Bowen  and  Inman) 1^,  and  rip  currents  (Bowen  and  Inman) 1^.  The 

lack  of  direct  observational  evidence  of  edge  wave  modes  on  beaches  (i.e.,  wave 

crests  pointing  offshore!)  naturally  raises  questions  as  to  how  and  if  these  modes 

can  be  excited.  The  most  important  mechanism  for  extensive  generation  opportunities 
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of  edge  waves  was  provided  by  Guza  and  Davis  who  demonstrated  theoretically 
that  Stokes ’  modes  could  be  excited  by  simple  wave  trains  normally  incident  onto 
a  perfectly  reflecting  beach  from  deep  water.  Basically,  the  nonlinear  coupling 
between  the  incident/ re fleeted  wave  of  frequency  W  and  background  edge  wave 
"noise"  of  subharmonic  frequency  u>  =  W/2  leads  to  a  resonant  interaction  and 
growth  of  the  edge  wave  noise.  Eventually  other  nonlinear  processes  develop 
to  limit  growth;  however,  at  steady  state  the  edge  wave  amplitude  exceeds  that 

of  the  incident/reflected  wave.  The  entire  evolution  of  the  edge  wave  has  been 

.  21 
investigated  theoretically  by  Guza  and  Bowen  (hereafter  referred  to  as  G-B) , 
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Minzoni  and  Whitman  (hereafter  referred  to  as  M-W) ,  and  Rockliff 

Herein  we  present  quantitative  experimental  measurements  on  edge  wave 
excitation  by  normally  incident  wave  trains.  Some  of  the  predictions  of  G-B  and 
M-W  are  tested  and  necessary  modifications  for  imperfectly  reflecting  beaches  are 
presented.  An  outline  of  the  presentation  is  as  follows.  In  Sect.  2  we  review 
the  theoretical  results  with  emphasis  on  the  prediction  which  can  be  explicitly 
tested  by  the  experiments.  In  Sect.  3  a  brief  description  of  the  experimental 


facilities  and  data  analysis  techniques  are  discussed.  A  comparison  of  measured 
and  predicted  data  appears  in  Sect.  4  followed  by  a  summary  of  the  major  results  in 
Sect.  5.  More  details  of  each  aspect  of  this  note  may  be  found  in  Lin  and  Hammack.^ 
2.  —  Review  of  the  Theories 

The  first  description  of  both  edge  wave  excitation  and  subsequent  evolution 
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was  developed  by  Guza  and  Bowen  using  the  shallow-water  equations;  hence,  their 

results  are  limited  to  small  beach  angles  (6  <<  1) .  With  clever  but  intuitive 

reasoning,  G-B  isolate  several  processes,  analyze  each  separately,  and  combine 

23 

linearly  to  yield  a  complete  evolution  model.  (More  recently,  Rockliff  has 

reproduced  some  of  the  G-B  results  relying  more  formally  on  the  governing 
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shallow-water  equations.)  Minzoni  and  Whitham  use  the  full  equations  and  provide 

the  most  formal  (mathematically)  description  of  the  edge  wave  evolution  process; 

their  results  remain  valid  for  arbitrary  beach  slopes  B  £  ir/2.  The  greater 

formality  of  the  M-W  formulation  permits  more  justifiable  ordering  of  different 

processes.  Hence,  the  outline  below  of  edge  wave  evolution  is  based  primarily 

on  their  formalism.  Following  the  qualitative  description,  a  quantitative  listing 

of  relevant  parameters  predicted  by  both  M-W  and  G-B  will  be  presented. 

The  classical  water-wave  problem  for  the  inviscid,  irrotational,  barotropic 

motions  of  an  incompressible  ocean  in  the  wedge-shaped  region  of  figure  11a  may 

be  formulated  in  terms  of  a  complex  velocity  potential  Initially,  we  have  a 

linear  wave  train  normally  incident  onto  the  beach  from  deep  water;  the  incident 

wave  amplitude  far  offshore  is  a^.  We  also  take  the  incident  wave  to  be  perfectly 

reflected  from  the  beach.  (The  reflection  coefficient  is  R  =  a  /a.  =1  where  a 

r  l  r 

is  reflected  wave  amplitude.)  Hence,  the  forcing  for  the  onshore  edge  wave  noise 
is  a  standing  wave  mode  which  is  a  member  of  the  continuous  spectrum  of  normal  modes 
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discussed  in  Sect.  1.  To  simplify  (somewhat)  the  standing  wave  description, 

it  is  further  assumed  that  the  beach  slopes  are  a  member  of  the  denumerably 

infinite  set  8  =  ir/2N  with  N  =  1,2,  ...  .  The  complex  potential  for  the 

standing  wave  may  be  written  as  $  '  exp(iwt) ,  and  we  note  that  ($  |  ~  e  <<  1. 

In  addition  to  the  primary  standing  wave,  we  assume  three-dimensional  perturbations 

by  Stokes  edge  wave  with  complex  potential  <j>e  *  exp(iwt)  where  |  <J>  |  <<  e  <<  1 

initially.  On  a  very  short  time  scale  (stage  1) ,  nonlinear  interactions  between 

the  primary  wave  and  the  edge  wave  noise  are  negligible.  However,  on  a  longer 

time  scale  (stage  2) ,  quadratic  interactions  arise  as  a  consequence  of  the  nonlinear 

boundary  condition  at  the  free  surface  and  become  significant.  Assuming  the  edge 

wave  noise  that  satisfies  the  resonance  condition  w  =  W/2  to  dominate,  quadratic 

interactions  <J>g  <J>e*  (where  *  denotes  complex  conjugate)  give  rise  to  terms  of 

the  form  exp(iWt/2)  and  contribute  to  the  growth  of  the  subharmonic  edge  wave 

noise.  This  interaction  between  the  linear  standing  wave  and  linear  edge  wave 

noise  produces  exponential  growth  of  the  form  a  ~  exp(yt)  where  a^  is  the  (real) 

edge  wave  amplitude  at  any  crest  location  (i.e.,  a  =  a(x,y,t) ,  and  we  take 

a(x,0,t)  =  A(x,t)  as  the  shoreline  amplitude)  and  y  is  the  initial  growth  rate. 

The  edge  wave  noise  grows  until  a  later  time  scale  (stage  3)  is  reached  where 

its  finite  amplitude  leads  to  the  development  of  processes  which  limit  further 

growth.  First,  quadratic  self -interactions  of  the  edge  wave  with  the  form 

<|>  4>  exp(iWt)  contribute  to  the  offshore  standing  wave;  this  is  termed 

6  6 

radiation  by  G-B  and  corresponds  to  the  fact  that  nonlinear  (even  second  order) 
standing  edge  waves  leak  energy  at  frequency  2u>  to  deep  water.  Since  1 4>s  |  ~  c 
for  the  disturbed  standing  wave  in  stage  3,  radiation  feedback  suggests  that 
a  steady  state  may  be  reached  when  |<J>e|  '  e*5 ,  i.e.,  the  edge  wave  becomes  larger 
locally  than  the  forcing  wave.  Second,  cubic  self-interaction  terms  of  the  form 


t>  <j>  d>  *  exp(iWt/2)  contribute  further  to  edge  wave  nonlinearity  by  modifying 
e  e  e 

its  dispersion  relation.  Hence,  the  system  is  retuned  and  the  resonant  frequency 
is  shifted.  If  the  forcing  was  perfectly  resonant  initially,  it  will  now  be  off 
resonance  and  an  effective  reduction  of  the  on-resonance  growth  rate  y  will  occur. 
(Of  course,  if  the  initial  forcing  is  slightly  off-resonance,  a  similar  modification 
in  y  is  required.)  When  combined,  all  of  these  processes  lead  to  an  evolution 
for  (inviscid)  edge  waves  of  the  form 

dA  2  2  **  3 

(41)  —  =  (y  -  k  )  A  -  y  A3 

where  y  is  the  initial  on-resonance  growth  rate,  k  is  a  measure  of  the  reduction 
in  forcing  efficiency  due  to  nonlinear  retuning,  and  y  is  a  feedback  coefficient 
due  to  both  radiation  and  retuning.  (Note  that  in  (41)  we  tacitly  assume  that 
the  edge  wave  phase  remains  constant  through  stage  3.  Theoretically,  this 
assumption  is  invalid;  however,  the  experimental  data  support  its  applicability.) 

At  the  risk  of  misquoting  G-B  and  M-W  and  possibly  introducing  numerical 
errors,  we  now  list  quantitative  expressions  distilled  from  their  studies  for 
the  parameters  appearing  in  (41) .  The  initial  (on- resonance)  growth  rate  predicted 
for  the  edge  waves  is 

0.0424  a.  W3 

(42a)  G-B;  y  =  - - — ' 

g/T  8 

0.0426  a.  W3 

(42b)  M-W  y  =  - 1 - ' 

g/S~  tan  B 

where  the  M-W  result  is  accurate  to  within  3%  for  B  ir/4.  These  predictions 
are  in  agreement  in  the  shallow-water  limit  (6  -*•  0)  .  Modification  of  the 
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initial  growth  rate  by  nonlinear  retuning  in  steady  state  may  be  conveniently 

written  to  the  same  order  of  small  parameter  £  =  ot  k  (where  a  is  the 

0  0  0 

steady-state  run-up  amplitude  of  the  edge  wave)  as 

....  „  (2)  (0) 

(43)  k  =  w  -  w 


with  <i/2^  and  representing  the  nonlinear  and  linear  natural  frequencies 

(2) 

of  the  edge  wave,  respectively.  The  shallow-water  limits  for  w  given 
explicitly  by  G-B  and  inferred  from  the  presentation  of  M-W  are 


(44a) 

G-B: 

(2) 

0) 

( 0 ) 

=  0) 

[i  + 

0.055  e  2] 

0 

(44b) 

M-W: 

(2) 

U) 

(0) 

=  .  0) 

[i  + 

0.012  £  2] ; 
0 

=  (g  k  B)*5  in  both  cases.  The  different  results  of  (44a)  and  (44b)  do 
not  appear  reconciliable.  Predictions  for  the  feedback  coefficient  y  in  the 
shallow-water  limit  also  differ  at  the  same  order  of  approximation  according  to 

W5 

(45a)  G-B:  y  =  0.001795  - -  ■ 

g  e 

(45b)  N-W:  u  =  0.002323 

g  s 


Finally,  we  note  that  the  steady-state  amplitude  of  the  edge  wave  at  the 
shoreline  is  found  from  (41)  to  be 

(46a)  A 
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2,4 

K  ) 
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3.  —  Experimental  Facilities  and  Procedures 

A  set  of  five  experiments  are  conducted  to  study  the  excitation  of 
standing  edge  waves  on  a  beach  by  waves  normally  incident  from  offshore. 
Experiments  are  performed  in  a  laboratory  basin  1.83  m  wide,  4.0  m  long 
and  30  cm  deep  as  sketched  in  figure  12.  The  tank  is  equipped  with  a 
mechanical  wave  generator  and  a  smooth  beach  whose  slope  can  be  varied.  An 
actual  beach  width  of  b  =  1.60  m  was  used  in  the  experiments  by  incorporating 
internal  sidewalls  over  the  beach  section  as  indicated  in  figure  12.  This 
enables  the  offshore  tank  section  of  uniform  depth  to  be  surrounded  by 
energy-absorbing  material  to  minimize  unwanted  reflections  and  simulate 
open-ocean  conditions.  Extensive  precautions  were  taken  to  insure  that  the 
shoreline  and  wavemaker  were  parallel  and  to  prevent  beach  "pumping"  during 
actual  experiments. 

An  offshore  water  depth  of  h  =  25.4  cm  was  used  in  all  five  experiments 
while  the  beach  slope  was  varied  from  6  =  15°  (N  =  6)  to  6  =  22.5°  (N  =  4)  ? 
one  experiment,  run  4,  was  conducted  at  6  =  20°  which  does  not  correspond  to 
a  slope  of  w/2N  for  any  integer  N.  It  is  immediately  apparent  that  the 
offshore  uniform  depth  region  does  not  conform  to  the  (unrealistic)  mathematical 
model  where  the  water  depth  increases  linearly  offshor  .  v  :ver,  the  offshore 
distance  over  the  sloping  beach  section  always  exceeded  three  e-folding  distances 
(>  95%  decay)  for  the  edge  wave  modes  excited  (assuming  Stokes'  modes).  Based 
on  Minzoni,25  it  is  not  expected  that  the  edge  wave  dispersion  relation  is 
affected  by  the  uniform  depth  offshore  in  any  of  the  experiments. 

In  all  experiments  wave  amplitudes  are  measured  at  three  locations;  one 
gage  onshore  over  the  sloping  beach  and  two  gages  offshore  in  the  uniform 
depth  region.  The  onshore  gage  is  always  positioned  at  a  longshore  location 
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corresponding  to  an  edge  wave  antinode.  The  reality  of  edge  waves  very 
near  the  shoreline  requires  that  they  deviate  from  the  profile  predicted 
by  (38) ;  hence,  the  onshore  gage  is  located  approximately  one  e-folding 
distance  away  from  the  shoreline  and  always  outside  of  the  surf  zone. 

The  procedures  adopted  for  each  of  the  five  experimental  runs  are 
as  follows.  First,  the  frequency  w  of  the  desired  edge  wave  is  calculated 
based  on  the  beach  width  b  =  1.6  m;  the  offshore  wavemaker  is  then  adjusted 
to  a  frequency  W  =  2w.  The  incident  and  reflected  wave  amplitudes  (a^  and  a^) 
of  the  partial  standing  wave  systems  are  then  determined  from  simultaneous 
measurements  at  the  offshore  gages  in  the  absence  of  edge  wave  excitation.  Edge 
wave  excitation  is  suppressed  by  inserting  a  thin  plate  perpendicular  to  the 
shoreline  which  penetrates  the  surf  zone;  the  plate  introduces  boundary 
conditions  which  destroy  the  resonance  condition  necessary  for  rapid  edge  wave 
growth.  The  thin  plate  is  then  removed  permitting  the  evolution  of  the  edge 
wave.  Once  the  edge  wave  attains  a  steady  state  the  wavemaker  is  stopped  and 
the  (unforced)  decay  of  the  edge  wave  over  the  sloping  beach  is  measured.  Hence, 
each  experiment  produces  a  time  series  of  the  water  surface  elevation  onshore 
and  offshore  containing  stage  1,  stage  2,  stage  3,  and  additional  stage  4  where 
edge  wave  forcing  is  terminated  and  viscous  damping  forces  dominate. 

4.  —  Comparison  of  Experiment  and  Theory 

Table  1  summarizes  the  measured  data  for  the  partial  standing  waves 
generated  by  the  wavemaker/beach  system  in  the  absence  of  edge  waves.  The 
incident  and  reflected  wave  amplitudes  are  presented  along  with  the  respective 
reflection  coefficients  which  range  from  3  =  0.134  to  R  =  0.500.  (Recall  that 


theoretical  results  assume  perfect  reflection  with  R  =  1.)  It  should  be 

noted  that  wave  breaking  near  the  shoreline  was  observed  in  each  experiment 

and  was  especially  strong  in  run  1. 

A  typical  time  series  (run  2)  taken  onshore  and  offshore  during  edge 

wave  excitation,  evolution,  and  eventual  damping  is  shown  in  figure  13.  The 

onshore  record  clearly  shows  the  effect  of  edge  wave  growth  as  it  alternates 

between  construction  and  destructive  interference  with  successive  c rests  of 

the  incident/reflected  wave  system.  (This  behavior  is  a  direct  consequence  of 

subharmonic  excitation  with  W  =  2a>.)  A  periodogram  computed  using  the  Fast 

Fourier  Transform  of  the  onshore  wave  record  for  run  2  is  shown  in  figure  14. 

The  first  peak  in  the  spectrum  corresponds  to  the  excited  edge  wave  mode  while 

the  second  peak  corresponds  to  the  standing  wave  generated  from  offshore. 

The  equally  spaced  peaks  at  higher  frequencies  represent  the  superharmonics  of 

both  the  edge  wave  and  standing  wave.  The  centering  of  wave  content  about 

well-separated  and  narrow  bands  with  identifiable  dominant  frequencies  (w,  2w,  3oj, 

such  as  that  exhibited  in  figure  14  permits  the  use  of  complex  demodulation 

techniques  to  view  the  real  time  evolution  of  both  the  amplitude  and  phase  of 
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each  dominant  wave  component  (see  Bloomfield  ) .  Using  the  periodogram,  an 

initial  guess,  say  w  ,  for  the  dominant  frequency  in  a  band  of  interest  is 

determined.  The  time  series  is  then  multiplied  by  exp(-iw  t)  to  shift  the 

a 

desired  wave  content  to  zero  frequency.  This  demodulated  signal  is  then 
low-pass  filtered  in  time  to  yield  the  instantaneous  phase  and  amplitude  of  the 
signal  component  with  frequency  cu  .  If  the  estimated  frequency  w  is  in  error, 

g  g 

a  linear  change  in  instantaneous  phase  with  time  will  be  observed;  the  slope  of 
the  linear  change  represents  the  error  in  In  this  manner  the  dominant 


frequency  of  each  narrow  band  in  the  periodogram  may  be  determined  very 

accurately  as  well  as  its  instantaneous  amplitude  and  phase.  (More  details 

on  the  application  of  complex  demodulation  technique  including  a  discussion 
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of  the  low-pass  filter  properties  may  be  found  in  Lin  and  Hammack.  )  The 

edge  wave  frequencies  in  each  experiment  were  determined  using  this  technique 

and  are  summarized  in  table  2.  Figures  15  and  16  show  the  instantaneous 

amplitude  and  phase,  respectively,  of  the  edge  wave  harmonic  at  the  onshore 

gage  in  run  2.  Both  parameters  oscillate  rather  wildly  during  the  initial 

time  t  <  25  secs;  this  behavior  is  a  characteristic  result  when  the  signal-to-noise 

ratio  is  small.  (Note  that  during  this  period  this  instantaneous  amplitude 

of  the  edge  wave  is  less  than  0.07  cmi)  The  edge  wave  amplitude  in  figure  15 

then  begins  to  grown  in  an  exponential  manner  while  the  edge  wave  phase  in 

figure  16  becomes  constant.  As  time  continues  to  increase,  edge  wave  growth 

slows  and  a  steady  state  is  achieved  for  t  >  70  secs;  the  gage-site  amplitude 

of  the  edge  wave  at  steady-state  is  aQ  =  0.754  cm.  The  measured  growth  rates 

Y  and  steady-state  amplitudes  a  are  shown  in  table  2.  The  inferred  shoreline 
in  o 

amplitude  Aq  based  on  the  measured  a^  at  gage-site  is  easily  calculated;  the 
ratio  of  A^/a^  ranges  from  3-6  in  all  experiments  in  agreement  with  the 
theoretical  prediction  that  steady-state  amplitudes  of  the  edge  wave  should 
exceed  those  of  the  offshore  standing  wave.  It  is  important  to  note  that 
measured  growth  rates  are  significantly  influenced  by  viscous  and  turbulent 
damping  forces  which  must  be  clarified  before  a  legitimate  comparison  with 
inviscid  and  irrotational  theoretical  models.  The  viscous  damping  rate  is 
easily  measured  in  the  laboratory  model  simply  by  turning  off  the  wavemaker 
and  monitoring  the  damping  of  the  edge  wave  harmonic  with  time.  Typical 
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results  are  shown  in  figure  17  which  corresponds  to  run  2.  As  expected, 

exponential  decay  occurs  when  the  edge  wave  becomes  sufficiently  small  (and 

linear  theory  is  applicable) .  The  measured  decay  rates  6  in  the  exponential 

stage  for  each  experiment  are  presented  in  table  2.  As  noted  by  Guza  and 
27 

Bowen,  the  damping  of  edge  waves  in  the  presence  of  breaking  incident  waves 

is  likely  to  be  dominated  by  turbulent  exchange  mechanisms  rather  than  (laminar) 

viscosity.  Viscous  effects  are  enhanced  by  laboratory  model  scales  and  should 

play  an  important  role;  however,  based  on  the  results  to  follow,  turbulent 

exchange  mechanisms  also  appear  to  be  important  at  laboratory  scales.  Ignoring 

the  effects  on  laboratory  damping  of  edge  waves  due  to  incident  wave  breaking, 

a  measured  "inviscid"  estimate  of  the  initial  growth  rate  may  be  calculated 

by  y  '  =  y  +  6 .  Results  for  y  '  are  tabulated  in  table  2  along  with  the 

m  ip  m 

theoretical  predictions  for  the  growth  rate  according  to  42a.  It  is  evident 

that  large  discrepancies  still  exist;  however,  it  should  also  be  remembered  that 

the  theoretical  predictions  assume  a  perfectly  reflecting  beach.  To  examine  the 

potential  effect  of  the  imperfect  reflection  in  the  experiments  we  have  plotted 

the  ratio  of  y  '  /Y  versus  the  reflection  coefficient  R  in  figure  18.  Note  that 

excellent  correlation  exists  (data  lie  on  45°  line)  suggesting  that  the  growth 

rates  should  be  calculated  using  the  reflected  wave  amplitude  a^_  instead  of  the 

incident  wave  amplitude  a^ ,  i.e.,  we  should  replace  y  by  yR-  Theoretical  estimates 

of  yR  are  listed  in  table  2.  It  is  important  to  emphasize  that  this  behavior 

of  the  measured  data  clearly  demonstrates  that  it  is  the  reflected  wave  component 

of  the  offshore  standing  wave  which  drives  the  edge  waves.  It  further  supports 

17 

the  observation  (Guza  and  Inman  )  that  edge  waves  do  not  occur  on  dissipative 
beaches  which  produce  little  reflection. 
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In  order  to  compute  the  theoretical  amplitude  of  the  edge  waves  in 
steady  state,  the  results  above  suggest  that  we  should  modify  the  theoretical 
growth  rate  in  (46)  to  the  form 

2  ,*r 

R(y  -  o 

V 

The  theoretical  amplitude  at  a  specific  gage  site  then  becomes 

(47)  a  =  A  cos  k  x  exp(-k  y  cos  B) ; 

oo  m  m 

results  of  computations  for  the  predictions  of  both  G-B  and  M-W  are  shown 
in  table  2.  Except  for  run  4,  theoretical  results  exceed  measured  data  by 
factors  of  1.18  to  1.91.  This  behavior  pattern  is  expected  since  the  turbulent 
damping  of  the  edge  waves  in  the  presence  of  breaking  incident  waves  has  been 
neglected.  Also  consistent  with  this  hypothesis  is  the  fact  that  agreement 
between  experiment  and  prediction  is  best  in  run  5  where  the  least  wave  breaking 
was  observed.  Run  4  is  anomolous  as  the  measured  amplitude  exceeds  theoretical 
prediction.  (Interestingly,  absolute  agreement  between  measured  and  predicted 
data  is  best  for  run  4!)  There  is  also  the  following  evidence  that  run  4  is 
unusual  re.  ive  to  the  other  experiments.  Note  from  table  1  that  the  offshore 
itiniina  wave  changes  little  between  runs  3  and  4  while  the  steady-state  amplitude 
•  •  wave  in  table  2  doubles!  In  other  words,  a  major  difference  in 

■  ■  -:r-  even  thouah  the  forcing  remains  approximately  constant. 
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ween  runs  3  and  4  is  the  slightly  greater 
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beach  slope  (2°)  in  run  4  which  no  longer  corresponds  to  one  of  the  set 
it/2N  required  by  the  theoretical  models.  Based  on  this  limited  datum,  it 
does  appear  that  the  predictions  for  steady-state  response  on  beach  with 
slopes  satisfying  6  =  t/2N  cannot  be  extrapolated  to  intermediate  slopes. 
(However,  extrapolation  for  the  initial  growth  rates  does  appear  permissible 
with  the  modifications  due  to  partial  reflection  described  earlier.) 

5.  —  Conclusions 

Based  on  the  experiments  described  herein,  the  following  major  conslusions 
may  be  stated  regarding  edge  wave  excitation  on  beac.nes  by  normally  incident 
waves  from  offshore. 

i)  The  theoretical  growth  rates  calculated  for  perfectly  reflecting  beaches 
may  be  generalized  to  beaches  with  imperfect  reflection  simply  by 
reducing  the  growth  rates  in  direct  proportion  to  the  reflection 
coefficient  of  the  beach.  This  result  suggests  that  the  reflected  wave 
is  indeed  the  driving  force  for  edge  wave  response  and  absolutely 
necessary  in  order  to  excite  these  modes. 

ii)  Growth  rates  modified  for  partial  reflection  effects  may  be  extrapolated 
to  beaches  with  slopes  not  satisfying  S  =  tt/2N. 

iii)  No  change  in  the  edge  wave  phase  occurs  once  the  exponential  growth 
stage  is  encountered.  Hence,  the  evolution  equation  for  the  real 
amplitude  of  the  edge  wave  is  applicable. 

iv)  There  is  limited  evidence  that  the  steady-state  amplitudes  of  edge 
waves  on  beaches  with  8  7*  ff/2N  are  significantly  larger  than  those 
for  nearby  beach  slopes  with  8  =  t/2N. 


i 


v)  Damping  effects  on  edge  waves  are  significantly  influenced  by  turbulent 
exchange  mechanisms  resulting  from  breaking  incident  waves,  even  on  the 
laboratory  scale. 

★  *  ★ 
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FIGURE  CAPTIONS 


Figure  1.  Definition  sketch  of  the  fluid  domain. 

Figure  2.  Schematic  drawing  of  the  wave  generator. 

2 

Figure  3.  Surface  wave  evolution:  h  =  5  cm,  b  =  61  cm,  V  =  30.5  cm  , 

N  =  3.  '  ■■  ,  measured  profiles;  #  ,  soliton  profiles 

computed  using  (7)  . 

2 

Figure  4.  Surface  wave  evolution:  h  =  5  cm,  b  =  30.5  cm,  V  =  30.5  cm  , 

N  =  3.  «—  ,  measured  profiles;  #  ,  soliton  profiles 

computed  using  (7) . 
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Figure  5.  Surface  wave  evolution:  h  =  5  cm,  b  =  61  cm,  V  =  30.5  cm  ,  N  =  0. 
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Figure  6.  Surface  wave  evolution:  h  =  5  cm,  b  =  30.5  cm,  V  =  -30.5  cm  ,  N  =  1. 

2 

Figure  7.  Surface  wave  evolution:  h  *  5  cm,  b  =  30.5  cm  ,  V  =  0,  N  =  2. 

—  hi  ■■  ,  measured  profiles;  •  ,  soliton  profiles  computed  using  (7) . 

Figure  8.  Internal  wave  evolution:  h^  =  45  cm,  h 2  =5  cm, 
b  =  61  cm,  A  =  0.05 

Figure  9.  Comparison  of  leading  wave  profiles  with  theoretical 
(internal)  KdV  soliton. 

Figure  10.  Evolution  of  a  narrow-banded  wave  packet.  -  »  ,  measured 

profiles;  #  ,  soliton  profile  computed  using  (12). 

Figure  11a)  Definition  sketch  of  fluid  domain  and  Stoke' s  edge-wave  mode. 

b)  Sample  offshore  profiles  of  higher  edge-wave  modes:  — —  , 

discrete  spectrum  (trapped)  mode;  — — ,  continuous  spectrum 
(untrapped)  mode. 
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Figure  12. 

Figure  13. 
Figure  14. 
Figure  15. 

Figure  16. 

Figure  17. 

Figure  18. 


Schematic  drawing  of  wave  basin  —  plan  view. 

Time-series  of  water  surface  elevations  onshore  and  offshore. 

Periodogram  of  onshore  time  series. 

Amplitude  evolution  of  subharmonic  edge  wave  for  run  2. 

Phase  evolution  of  subharmonic  edge  wave  in  run  2. 

Viscous  damping  of  subharroonic  edge  wave  in  run  2. 

Correlation  of  measured  and  theoretical  "inviscid"  growth 
rates  with  beach  reflection  coefficient. 
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